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1. INTRODUCTION 

The Ising-like nature of the universal features of the 
one-component fluids close to their vapor-liquid critical 
point is now well-estabHshed (author?) These uni- 
versal features can be estimated by using a renormaliza- 
tion approach (author?) [2| of the classical-to-critical 
crossover behavior (author?) for three-dimensional 
(3D) Ising-Hke systems with a symmetrical order pa- 
rameter density. Especially, the massive renormalization 



scheme (author? 



lispecially, 

?) i, a a 



applied to the $^^3 (n = 1) 



model {d and n are the dimensions of the space and order 
parameter density, respectively) , has been recently revis- 
ited (author?) [7[ to provide max and min crossover 
functions which include updated estimations (author?) 
[1| of the universal values for exponents and amplitude 
combinations. Subsequently, a related paper (author?) 

has provided a unique universal form of the mean 
crossover functions, which enabled asymptotic descrip- 
tion of the universal features valid within the Ising- 
like preasymptotic domain exactly, by eliminating the 
increasing uncertainties associated with the asymptotic 
error-bar propagation in the initial min and max estima- 
tions of the critical exponents. 

These crossover functions have been used to analyze 
for example, the singular behavior of the correlation 
length of seven pure fluids in their homogeneous do- 
main (author?) or the singular behavior of the 
squared capillary length of twenty pure fluids in their 
non-homogeneous domain (author?) For both 

cases, it was shown that the scale dilatation method pro- 
posed by one of us (author?) [12, HI, [11, [H, ^ gives 
a corresponding master singular behavior for the one- 
component fluid subclass which agrees with a description 
by an appropriate modiflcation of the theoretical func- 
tion (author?) [l3| within the Ising-like preasymptotic 
domain. As a matter of fact, such an exact behavior of 
each theoretical function, approximated by a two-term 
Wegner-like expansion (author?) [IS] close to the non- 



trivial flxed point, is an essential tool to provide better 
understanding of the basic relations between the "mea- 
sured" asymptotic amplitudes and the unknown (fluid- 
dependent) scale factors introduced by linear approxima- 
tions between bare flelds and physical flelds (author?) 
[l9| . From the analogy with the scale dilatation method 
which relates the master flelds to the physical flelds intro- 
ducing two well-deflned critical parameters of each fluid, 
it was then recently proposed (author?) Jl3| an unam- 
biguous modiflcation of the mean theoretical functions 
to represent the master crossover for the one-component 
fluid subclass (labeled {!/}) in conformity with the corre- 
sponding universal features of the Ising-like universality 
class (labeled $3 (1)). In the following, the theoretical 
functions (valid for any three-dimensional Ising-like sys- 
tem) estimated in Ref. (author?) Q are called mean 
crossover functions and the related reference (author?) 
Q is labeled I, while their modiflcation (only valid for 
the one-component fluids) proposed in Ref. (author?) 
[l7| are called master crossover functions and the related 
Ref. (author?) ^ is labeled II . 

The vaHdity range of the mean crossover functions is 
theoretically funded only within the Ising-like preasymp- 
totic domain (i.e. for t < Cp™^, where t is the thermal- 
like fleld along the critical isochore, i.e. for h = 0, 
where h is the magnetic-Hke fleld , see II). Especially, 
the uniqueness of the asymptotic scale factor, which 
acts as a crossover parameter along the critical isochore, 
has never been clearly demonstrated for the effective 
extended asymptotic domain where the mean crossover 
functions flt correctly the experimental results. The main 
goal of the present paper addresses to the explicit calcu- 
lation of this single crossover parameter using the mean 
crossover functions beyond the Ising-like preasymptotic 
domain. 

We consider in detail the singular behavior of the 
isothermal compressibility kt (At*) and the correlation 
length ^(Ar*) of xenon as a function of the reduced 
temperature distance At* = along the critical 
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isochore, in the homogeneous domain T > Tc, [T (Tc) 
is the temperature (critical temperature)]. Xenon is 
here selected as a standard "critical fluid" in the sense 
that the present work only uses quantities which orig- 
inate from theoretical arguments, now well-understood, 
to link three independent dimensionless amplitudes to 
the needed three scale factors (author?) [ioj- As a mat- 
ter of fact, the fit of susceptibility and correlation length 
data obtained by Giittinger and Cannell (author?) \2l\ 
from their precise turbidity measurements in the tem- 
perature range 0.028 X < T -Tc < 29K, can then be 
used as representative of an "ideal" result in order to 
check carefully the temperature range where the value 
(and its attached uncertainty) of the scale factor char- 
acteristic of the thermal field is effectively determined. 
This singular behavior was first analyzed (author?) \2^ . 
jointly with singular behaviors of the correlation length 
(author?) fill the heat capacity (author?) 
using precisely the initial estimations (author?) Q of 
the crossover functions from the massive renormalization 
scheme. In addition to the experimental critical tem- 
perature Tc{~ 289. 74K), pressure Pc (- 5.84 MPa), and 
density pc (— 1110 kgm~^), for the first time a minimal 
quantity made of three Ising-like non-universal parame- 
ters of xenon was introduced as a set of a single critical 
wavelength (noted go) and two dimensionless scale factors 
(noted ■& and ip) for the thermal-like (t) and magnetic- 
like (h) fields, respectively, with the analytical relations 
t — i?Ar* and h = i^pAfl vaHd for the asymptotic Hmits 
t -> and h ^ (author?) [19, 2g]. A/i is the stan- 
dard notation of the dimensionless ordering field in the 
fiuid case (author?) [l^l, which is related to the chemi- 
cal potential difference to the critical chemical potential 
[see below Eq. (flGl) ]. In the following, when we refer 
to this initial result obtained by fitting the experimen- 
tal data with the (max) crossover functions estimated for 
the sixth-order, we will use the label MR6niax (author?) 

aa. 

Now, using the updated mean crossover functions given 
in I, the main objective is to replace the fitting adjust- 
ment of the fiuid-dependent parameters go, and ipp, 
by their exact asymptotical values defined in II, only 
introducing three master numbers (noted 9^^^^, L'f^-^^, 
and ^''fi^}) which are characteristics of the {l/}-subclass 
(see II). In that scheme, go-, i?-, and ^pp-vahles origi- 
nate from calculations of the three fluid-dependent ampli- 
tudes , r~^, and a+ (in standard notations (author?) 
[23. l25|l. using the scale dilatation method where each 
fluid is characterized by only four well-defined critical 
point coordinates. Such results account for the master 
singular behavior of the {l/}-subclass observed within 
the Ising-like preasymtotic domain bounded by a known 

limit (defined as >Cp^]5 = ^ijj in II) , in conformity with 
the two-scale factor-universality of the $3 (l)-universality 
class. As a main result explicited below [see Eqs. l(35|) 
and l(40|)]. we unambiguously probe that the tempera- 
ture range Ar^j^^ < At* < Ar^^^ covered by the Giit- 
tinger and Cannell's experiments is beyond the Ising-like 



preasymptotic domain, i.e. Ar*^jj^ (~ 10~ ) > . 

Moreover, after the initial analysis of Ref. (author?) 
[2^ . the Giittinger and Cannell's xenon data were also 
used (author?) 26. 27, 28, 29, 3^ in support for discus- 




sion of several theoretical (author?) 1311. 1 32 
and phenomenological (author?) [3g.l37l.l38 

proaches of the crossover phenomena. As an essential 
common result, a single temperature-like crossover pa- 
rameter seems appropriate for a complete characteriza- 
tion of the classical-to-critical crossover in pure fluids, 
whatever the selected crossover theory or the phenomeno- 
logical approach (for comparative analyses, see for exam- 
ple Ref. (author?) and references therein) . However, 
the theoretical understanding in terms of a Ising-Hke crit- 
ical crossover characterized by the asymptotic value of 
the scale factor 1} appears limited for two main reasons: i) 
the fltted data are not within the Ising-like preasymptotic 
domain; ii) any crossover function is only well-deflned 
within the Ising-like preasymptotic domain t < Cp™^ 
(see I and II). Especially for the susceptibility case, such 
a result cannot be used in comparative analyses which 
debate on the "correct" shape of the temperature depen- 
dence of the effective exponent 7 defined as universal 
functions of a sin gle dimensionless variable (author?) 

[H, [13, [H, M, MM, M, 

Our present attention is then mainly focused on the 
determination of a thermal- like scale factor, noted 
calculated from the application of the mean crossover 
function beyond the preasymptotic domain. Our analytic 
procedure developed hereafter aims to retrieve the Ising- 
like universal features only using well-defined energy and 
length units, and a well-controlled number (three) of well- 
defined dimensionless scale factors given in a form of a 
set §si? as defined in II [see also Eqs. (|25|) and (|26|) be- 
low]. We strictly avoid adjusting the system-dependent 
parameters by a minimization of fitting errors to demon- 
strate that dc = Alternatively, using both estima- 
tions of the effective exponent and amplitude values (au- 
thor?) we calculate the local value of dc (Ar*) to 
verify the uniqueness of the -de value and its identity with 
the d value when Ar* 0. Accordingly, our method 
offers the great advantage to directly re-use previous re- 
sults obtained by fitting experimental data with an ef- 
fective power law valid in a limited temperature range 
(author?) [13, [12, [i3| • Therefore our attention is also 
focussed on the corresponding analyses of the isothermal 
compressibility data obtained from pVT measurements 
(author?) [3, [H, [13, [13, [Hi covering the temperature 
range 0.1 K < T — Tc < 283.41 K, light scattering mea- 
surements (author?) [2l|, [il, [5^, [sil, [13, [13| covering the 
temperature range 2.6 mK < T — Tc < 29 K, and inter- 
ferometry measurements (author?) [53,[1E,[H^ covering 
the temperature range —0 < T — Tc < 29 if , extending 
then significantly the temperature range investigated by 
the present study. We simultaneously provide the effec- 
tive reduced temperature range of the extended asymp- 
totic domain, bounded by a limit noted >Cead' with 
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(see II and below) , where isother- 



-^EAD ^ -^PAD ~ i5(Xe) 

mal compressibility of xenon is accurately described by 
the mean crossover function for susceptibility using the 
single scale factor (Xe) for the thermal-like field. Such 
a result gives the first complete comparison between ex- 
perimental results and crossover theories, de facto with- 
out any adjustable parameters in a relative temperature 
range covering more than four decades, when the mean 
crossover functions are appropriately modified to account 
for master behaviors of the one-component fiuid subclass. 

The paper is organized as follows. 

In Section 2, we recall useful notations and definitions 
needed to use the mean crossover functions for the cor- 
relation length and the isothermal susceptibility in the 
homogeneous phase of a one-component fiuid. In Section 

3, the characterization of the Ising-Hke preasymptotic 
domain of xenon is analyzed only using the isothermal 
compressibility fitting result obtained by Giittinger and 
Cannell, from their turbidity measurements performed 
outside the Ising-like preasymptotic domain. In Section 

4, after the introduction of the three-parameter charac- 
terization when the extension of the crossover domain 
remains undefined, we demonstrate the great advantage 
of the mean crossover functions to provide unambigu- 
ous scaling determination of dc beyond the Ising-like 
preasymptotic domain. 

In Appendix A, we give the basic estimation of the 
needed amplitudes from application of the scale dilata- 
tion method to xenon. All the other measurement meth- 
ods of the isothermal compressibility are then considered, 
with a special attention to the data obtained from pVT 
and interferometry measurements which infer a practical 
three point calibration for their relative comparison to 
the Giittinger and Cannell's data. We finally compare in 
a detailed manner all the results to the ones calculated 
with the mean crossover function for susceptibility, then 
justifying a posteriori the exact values of the three scale 
factors 5o, and V'p for xenon which are used in this 
paper. 



2. MEAN CROSSOVER FUNCTIONS FOR 
CORRELATION LENGTH AND 
SUSCEPTIBILITY 

2.1. Definitions and notations 



where t (> 0) is the thermal field like variable, h (= 0) 
is the magnetic field Hke variable. D {t) is a univer- 
sal crossover function for the confiuent exponents A and 
Amf given by 



D{t) = 



AMF>5'2%/t + A 



(3) 



All the universal exponents 7, A, Amf, theoretical 
amplitudes Z^, Z+, and theoretical parameters S2, ^^i, 
Y^^, A"^j, Y^-, are defined in I. They are also reported 
in Table [J for use in the following numerical estimations. 

More generally, the crossover functions estimated in 
I are only well-defined for the critical line that links the 
Gaussian fixed point and the non-trivial fixed point. The 
effects due to the second-order (and higher) analytical 
contributions and the ones due to the confiuent correc- 
tions to scaling linked to critical exponents A2, A3, 
have been discarded. As a direct consequence, these func- 
tions account exactly for the Ising-like universal features 
only estimated in the t-range very close to the non-trivial 
fixed point, which corresponds the Ising-Hke preasymp- 
totic domain t < ^Cp^D defined in I [see also below Eq. 
©]■ 

Within the Ising-like preasymptotic domain, the mean 
crossover functions can be approximated by their re- 
stricted (two-term) Wegner like expansion (author?) 
Then Eqs. ^ and ^ can be replaced by: 



t.PAD,th 



it) 



XPAD.th (i) = (Zj) 'i-^[l + Zi'+t^] 



(4) 



(5) 



where the amplitudes Z^'"*" and Z^^+ of the first-order 
term due to the lowest confiuent corrections to scaling 
are given by 



4-" = -JlK.y. 



+ 



(6) 



(7) 



The dimensionless mean crossover functions 
Fg {t, h = 0) = for the inverse correlation length, 

and {t, h = 0) = ^^{t) ^'-'^ inverse susceptibility, 
in the homogeneous phase, read as follows (see I), 

3 ^+ 



i=l 
3 

[xth (t)]"'=zji-^n(i+^xV 



(*) 



(1) 

(2) 



(see Table [D. They are related by the universal ratio 
(author?) 



Z^'+ 

-4-T = 0.67919 



(8) 



The validity domain of Eqs. (0]) to ([8]) corresponds to 



t < >CpAD where 



i=l 



Rising / 0-033 

PAD - \ 



= 1.9 X 10" 



(9) 



4 



a 


exponent 




& 


i 






b 


exponent 


Z+ 




i 






V 


0.6303875 


2.121008 


22.9007 


1 


40.0606 


-0.098968 


7 


1.2395935 


3.709601 


22.9007 


1 


29.1778 


-0.178403 


A 


0.50189 






2 


11.9321 


-0.15391 


A 


0.50189 






2 


11.7625 


-0.282241 


Amf 


0.5 






3 


1.90235 


-0.00789505 


Amf 


0.5 






3 


2.05948 


-0.0185424 














5.81623 












Zi'+ = 


8.56347 



Table I: Values of the universal exponents and constants of Eqs. {TJ to ([Zl), for (a) the dimensionless correlation length and (b) 
the dimensionless susceptibility, in the homogeneous domain (see Ref. (author?) Q). 



with S2 = 22.9007 (see I and Table [H 

Finally, any restricted two-term Wegner-like expan- 
sions can be estimated only using the following set of 
three theoretical amplitudes (author?) (20| . 

[t < C^ll:h = 0) = (Z+)"' ; (ZJ)-^} 

(10) 

where the subscript A labels for the amplitude nature of 
these characteristic parameters. 

Beyond the Ising-like preasymptotic domain, i.e. t > 
£p™p, we also recall that the theoretical parameter S2 
acts as a convenient sensor to estimate some orders of 
magnitude of t which are convenient for the analysis of 
the crossover. For example, the theoretical crossover 
ternperature Ia, which corresponds to the value D (1^) = 
""^ of the universal confluent function of Eq. ^ , is 
defined by 

tA^ (^^^ ^ 1.9 X 10-3 (11) 

For the confluent corrections, tA characterizes the uni- 
versal crossover exchange between predominant Ising-like 
nature close to the non-trivial fixed point [t ^ ^a], to 
predominant mean field-like nature close to the Gaussian 
fixed point [t ^ Ia]- Indeed, tA gives an estimation of the 
order of magnitude of the t„„ 1 -value where each effective 

exponent ep^e,th (t) — — ^^'/Inf/)*"'^ (author?) [4l| crosses 

its "mean" crossover value ep i (^t^^ 1^ — ^'"^^^•^^ (see 

Figure 4 in I). In such a "critical-to-classical" crossing 
range in the homogeneous domain [t > 0], 4h — 20 — 30 is 
a typical (dimensionless) order of magnitude for the the- 
oretical correlation length. For example, 4h (^a) = 28.8, 

Ah {^v,^ ~ 22.2, and ^th (^7,5) ~ 18.8, anticipating the 
following discussion of the results reported in Table Ull 
where the values of Eqs. (U) and ^ are calculated for 
eight conditions of the effective exponents. Moreover, in- 
troducing the practical relations >CpAD = '^^A, with w = 
10~^ (see I), it is easy to separate the analysis of, either 
the Ising-like preasymptotic domain t < w x Ia, or the 
intermediate "Ising-Hke" crossover domain tut a < t < Ia- 
Especially considering the selected values of the effective 
exponent 7e,th (t) = given in line 3 of Ta- 

ble |ll]), we note that the conditions [tfi ) (column 1) 



and 7e [t^^ j (column 2) are obtained for t > Ia, while, 

obviously, the condition 7e (^a) (column 3) is obtained 
at t = tA, where 4h reach a value ~ 30. For 4h ^ 30, 
or t < t A 1 , (columns 4 and right) , we expect to observe 

the Ising-Hke universal features of the critical phenom- 
ena which are then characterized by the three system- 
dependent parameters only asymptotically well-defined 
when £th ^ 1916, or t < Cp'^^ (i.e. within the Ising- 
like preasymptotic domain). We note that the effective 
values J^n^pVT = 1.16665 (author?) [Ti|), 1.19 (au- 
thor?) [2i,[4a|, 1.211 (author?) [42,] obtained from pFT 
measurements, are precisely in this intermediate Ising- 
like range 

'^PAD <t <tA. So that, we have underlined 
these three latter conditions in columns 5 to 6, respec- 
tively, showing that the value of theoretical correlation 
length £th is on the range 50 - 220, i.e., 2 < j^^ff^ < 8. 
A special attention to the finite "Ising-like" temperature 
range covered by pVT experiments is given in Section 
4.2 and Appendix A3. Correlatively, within the Ising- 
like preasymptotic domain t < £p]J^]f, the theoretical 
values reported in lines 7 and 8 of Table |ll] show that 
the condition £th > 2000 (or j;^^ > 70) is satisfied. 
The expected temperature- like variation of 7e,th (t), i.e., 
7 — 7e,th (-Cp^Q^ < 0.006, is then significantly lower than 
the typical uncertainty on the asymptotic experimental 
value 7exp ^ 1.23 ± 0.02 (author?) (H. 

The eight indexations of t (or ^th) given by these eight 
conditions of Table [ill will be used to label the horizon- 
tal (upper) axes of the next figures with corresponding 
arrows and column numbers. 



2.2. Fluid-dependent parameters 

The introduction of the fiuid-dependent parameters 
to fit the experimental results using the mean crossover 
functions was detailed in § 3.1 of II. We recall here the 
main steps to introduce definitions and notations of the 
needed quantities, fixing the value A*^ = 1 to neglect 
the quantum effects at the microscopic length scale (au- 
thor?) for the xenon case. 

When t ^ at ft. = 0, t is analytically related to 

T — T 

At* = — ^ -> (12) 
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(a) \ label 


1 


2 


3 


4 


5 


6 


7 




Condition 


/ /MF 

- 2 


- 2 


- 2 


7e,pVT(author?) fl2] 


7e,pVT (author?) [42] 


7e,pVT (author?) [24] 


•^PAD 




Te.th 


1.1198 


1.12776 


1.14081 


1.16665 


1.19 


1.211 


1.23397 






0.561 


0.5652 


0.572 


0.5862 


0.5994 


0.612 


0.6266 






0.500768 


0.500834 


0.500945 


0.501176 


0.501396 


0.501601 


0.50183 




t 


4.06 X 10"^ 


3.056 X 10"^ 


1.907 X 10"^ 


7.033 X 10""' 


2.392 X 10""* 


6.22 X 10"^ 


1.907 X lO""^ 






18.76 


22.0 


28.8 


51.32 


97.29 


220.1 


1916 






0.73735 


0.70492 


0.65161 


0.54730 


0.45641 


0.37810 


0.29363 




(b) 




















0.193 


0.145 


0.0905 


0.0334 > £ilo 


0.01135 


0.00295 


0.905 X 10"'' 




« _ ttth 

ar n,{l/} 


0.730 


0.857 


1.12 


1.997 


3.786 


8.564 


74.56 




r+ 

e,th 


0.099601 


0.098193 


0.095456 


0.088586 


0.080842 


0.072626 


0.061629 





Table II: (a): Calculated theoretical values of the effective exponents 7e,th (line 3), fe.th (line 4), A^.th (line 5), the thermal-like 
field t (line 6), the correlation length i'th (line 7), and the effective amplitude Zjg (line 8) of the susceptibility, using Eqs. ^ 
and ©J for the typical conditions given in line 2, respectively; (b): Corresponding physical values of the reduced temperature 
distance Ar* = |- (line 10), the dimensionless correlation length ^ = jjrh" C'^^ ■'■■'■)' ^'^'^ effective amplitude T+^j^ (line 
12) of the dimensionless compressibility, using Eqs. iflSll . l[2T|l . and lf23|l. for the xenon case, with i9 = 0.21069, L^^-'} = 25.6988 
and ipp = 3.2507 10"" (see text). 



by the following linear approximation (author?) [l9| 

t = i?At* (13) 

The dimensionless scale factor i? is a fluid-dependent pa- 
rameter. The fluid critical temperature Tc, provides the 
energy unit 

(/3c)"' - fcsTe (14) 

which links the dimensionless free energies of the 
model and the selected fluid. Accordingly, the dimension- 
less form of the Hamiltonian of the $''-model results of 
the introduction of a flnite (but arbitrary) wave number 
Aq, the so-called cutoff parameter (see I), whose inverse 
is related to the unknown flnite short range of the mi- 
croscopic molecular interaction. A convenient method at 
d — i consists in replacing (unknown) Aq by go which 
is the adjustable critical coupHng constant of the 
term having correct wave number dimension. The in- 
verse (go) acts as the physical length unit to link the 
theoretical dimensionless correlation length (^th) and the 
physical correlation length (^exp) of each one-component 
fluid, through the fltting equation 

4h (t) = .9o6xp (At*) (15) 

The dimensionless ordering fleld for fluids is deflned as 
(author?) 

Afl^ ftp- jlp^c (16) 

where A/l is written using practical dimensionless chem- 
ical potential jl = (author?) ji^j. Hp (Mp.c) is the 
(critical) chemical potential per mass unit, p {pc) is the 
(critical) mass density, and p (pc) is the (critical) pres- 
sure. Correspondingly, the practical dimensionless form 



Ap of the order parameter density reads as follows 

(17) 

Pc 

using the practical dimensionless form p = j- of the mass 
density (author?) 

When h ^ and A/2 ^ at t = 0, ft. is related to Ap, 
as follows 

h = ippAp, (18) 

where ipp is the second (fluid-dependent) scale factor. 
However, the parameters Tc, go, i9, and ipp, do not pro- 
vide unequivocal link between theoretical and physical 
thermodynamics quantities. For example, the deflnition 
of Po introduces a second unit — ~ r energy i £ (gpg- 

Pc L mass J ^ ^ 

ciflc) energy which differs from {mp(3c)~^ by the factor 
Zc [ mp is the mass of the fluid particle (i.e. the molecular 
mass)]. Correlatively, the critical mass unit of the one- 
component fluid introduces a critical speciflc volume -j- 

which differs from the volume I^s^jl = (ac)'^ of the critical 
interaction cell [see below Eq. ifTO]) ]. The comparison be- 
tween the two volumes introduces the extensivity of the 
system through the amount ^ of matter fllling the vol- 
ume of the critical interaction cell (author?) |l3|. That 
provides alternative choice between two energy units and 
two length units originating from thermodynamics. Such 
noticeable differences in the system-dependent units of 
the dimensionless variables impose to have careful atten- 
tion when comparing the dimensionless thermodynamics 
potentials. Here, the selected length unit ac reads 



6 



and have physical meanning in terms of the range of in- 
termolecular interaction in fluids. Therefore, in our no- 
tations, the superscript star labels a dimensionless vari- 
able which uses (/3c) and Uc as energy and length units 
in one self-consistent procedure to made dimensionless 
all the thermodynamic variables normalized per particle 
(author?) not per mass unit. The subscript p 

recalls for practical order parameter density defined by 
Eq. (flTl) . and the related practical dimensionless vari- 
ables are decorated by a tilde. For example, the experi- 
mental isothermal susceptibility XT,p for fluids is deflned 

by Yt o = I I ^ z — i when the fluid or- 

j ,fj yd^p J rp [ energy X volume J 

der parameter (respectively, the fluid ordering fleld) is 
proportional to the (mass) density p ~ [ volume ] 
spectively, the chemical potential per mass unit fip ~ 
r energy 1 N j.gjg^^gj ^j^g isothermal compress- 

ibihty nr = i (|)^ ~ [^] by xr.p = p'nr. 
The practical fluid dimensionless variables are p = j^, 
p = while = PcKT (author?) Q. Thus we 

obtain xt = '«t O'^ly at p = 1 for At* > 0. Such 
a practical interrelation between the dimensionless vari- 
ables results from the implicit use of these two "thermo- 
dynamic" length units previously defined from two dis- 
tinct volumes v- ^ = ^ and (ac)'^. On the other hand, 
from the theoretical scheme applied to the $^^3 (n = 1)- 
model, after normalization of the free energies by ksT = 

fc^Tc, the length dimensions of the Hamiltonian quanti- 

2 - - 

ties ro - Tqc ^ (go) t, ~ (50) ' m, and ho ~ {go)^ h, 

lead to xo.th = (50)^ Xth ^ [i^S^] (^^^ ^ notations 
and definitions). We can then conclude that introducing 
the dimensionless scale factors i? and ipp through Eqs. 
(fT3|) and ifTS)) , provides a subtle critical combination (not 
discussed here) of the model units and the fluid units. 
The essential point is to guarantee the uniqueness of the 
energy unit and the length unit in the description of di- 
mensionless singular behaviors (author?) [23|. Our se- 
lected length unit ac [see Eq. (fT9| ] takes thermodynamic 
origin, while the wave number unit go [see Eq. I|15p ] has 
theoretical interest to flt the asymptotic singular diver- 



(At* 



Then, by exchanging 



gence of Cxp 
Eq. (fT5|) and the following dimensionless form 



£th(t)=L^i^>Cxp(Ar*), (20) 

we also introduce the supplementary scale factor 

L^i^> = acgo (21) 

as a dimensionless product between the two critical quan- 
tities (i.e., deflned for i = 0; ft, = and At* = 0; Ap = 0, 
respectively). 

After all, m ^ and Ap — > are related by the equa- 
tion 

m=(Li'fA'\i;,y'Ap, (22) 



where m is the theoretical magnetization-like order 
parameter. Considering the theoretical susceptibility 



Xth (t) = (^)f and the experimental isothermal suscep- 
tibility XT,exp(Ar*) = (^f^^^ for fluids, the second 
fltting equation is then obtained as follows 



Xth (0 - (l^'^^) '(V^p)"'xT,exp(Ar*) 



(23) 



with XT,exp (Ar* 



K. 



T,exp 



(Ar*) when Ap = 0. 



Finally, each one-component fluid is asymptotically 
characterized by the set (author?) |! 2o| 



Q4Ar* ^0) - [{(3c) ';ac;go;^;^p} 
which can be rewritten in an equivalent form 

Qc(Ar*^0) = {(/3e)"';ac;§5f} 



(24) 



(25) 



then introducing the following fluid set Ss p made of three 
(asymptotic) dimensionless scale factors 



(Ar*^0) = {79;L{i^>;^p} 



(26) 



The subscript SF recalls for the scale factor nature of 
the three fluid-dependent dimensionless numbers, while 
the condition Ar* indicates the asymptotic nature 
of the hypotheses needed by the renormalization scheme. 
Within the Ising-Hke preasymptotic domain, (Ar*) 
and ^ (Ar*) can be approximated by 



T.exp 



(Ar*) =r+(Ar*)"^ 1 + aj (Ar 



^exp (Ar*) _ ^+ 



= ^ (Ar* 
ac 



l + a+ (Ar*)' 



with (author?) [3| 



Z 



1,+ 



Zi 



0.67919 



(27) 



(28) 



(29) 



Hereabove, we have selected r+, (or ^+ = ^), and 
aj as independent amplitudes to characterize each one- 
component fluid. The "experimental " parameter set writ- 
ten as 



ic 



I (/3c) ^■,ac;SAj 



(30) 



is then Ising-like equivalent to the set Qc {At* 0) of 
Eq. l(25l) [here the subscript /^pXd recalls for the fluid- 
dependent temperature domain of validity]. Its dimen- 
sionless part Sa (where the subscript A recalls for the 
amplitude nature of the dimensionless numbers) is given 
as : 



(31) 



which is compared with Eq. l(26|) using the fltting equa- 
tions (|23l) and l(20|) [or (flSl) ]. The successive unequivocal 



7 



determinations of the scale factors, first hence L''^-''^ 
(or go), and finally Vp, are obtained from the following 
hierarchy of equations: 



(32) 



or go 



r+z+^^ 



(33) 



(34) 



Each one-component fluid characterized by Qc of 
Eq. l(25l). has the Ising-Hke universal features of the 
^d=3 ~ l)-model in the Ising-like preasymptotic do- 
main >Cp^j3 given by 



3. XENON CHARACTERIZATION WITHIN 
THE ISING-LIKE PREASYMPTOTIC DOMAIN. 

3.1. Hypothesized description of the Ising-like 
preasymptotic domain 

We hypothesize that the two terms of the asymptotic 
Wegner-Hke expansion of K^.exp i^'^*) and^e^ip (Ar*) are 
exactly known for critical xenon, provided by the appli- 
cation of the scale dilatation method given in Appendix 
A. In such a situation, all the needed information takes 
origin on four critical coordinates [see below Eq. l|A2p ] 
which localize the xenon critical point on the experimen- 
tal phase surface of equation ^{p,vp,T) = 0. Accord- 
ingly, the (dimensional) values of Eqs. (fT4l) and (flOl) are 
the following 



Tc = 289.733 K or (/JJ"^ = 4.0002 x IQ-^i J 
ac = 0.881508 nm 

leading to the dimensionless amplitude set 



(38) 



^ •'-PAD 



/rising ^ 



1.9 X 10-6 
t9 



(35) 



Equation (|35|) demonstrates that the knowledge of the 
temperature-like scale factor ■& defines the extension of 
the Ising-like preasymptotic domain of each selected 
fluid, then providing an essential tool for analyzing ex- 
perimental data. Correlatively, admitting a single d value 
whatever the property and the thermal field range, the 
(fiuid dependent) crossover temperature Ar^ is given by 



tA 



1.9 X 10" 



(36) 



and the crucial problem of how to define the tempera- 
ture range At* < CpXi) solved. Appropriate rewriting 
of Eqs. l(32|) to l(34|) . provides the following functional 
scaling form 



Sa^{Ssf) with At* < 



1.9 X 10" 



(37) 



where F (§sf) are universal functions. We note the "the- 
oretical" (i.e. originating only from the MR scheme) na- 
ture of the l.h.s. of Eq. (|37)) . 

In next Section 3, we analyze the status of this ex- 
pected three-scale-factor characterization for the xenon 
case, first, using the values inferred from the applica- 
tion of the scale dilatation method (author?) [iSL and 
second, comparing the results (author?) [26, _23| ob- 
tained from the massive renormalization (MR) scheme 
(present work), the minimal subtraction renormalization 
(MSR) scheme (author?) [29|, and the crossover para- 
metric model (CPM) (author?) [i^l, with the experi- 
mental results of Giittinger and Cannell. 



Sa 



with 



a+ = 1.23397 



e+ = 0.209111 {^^ = 0.184333 nm) } (39) 
r+ = 0.0578204 



0.83810. Obviously, the vaHda- 



tion of our hypotheses and the justification of Eq. l(39|) 
require the detailed analysis of the isothermal compress- 
ibility data of xenon given in Appendix A. However, we 
note that the essential aspects for the following presen- 
tation are the Ising-like nature (since the fiuid character- 
ization originates from its critical point coordinates) and 
quantity (three) of the dimensionless amplitudes, while 
the quoted precision of their numerical values is of sec- 
ondary importance. Using Eqs. l(32|) to (|34l) . the three 
dimensionless scale factors for xenon are: 

d = 0.021069 
SsF = { Lii-'^} = 25.6936 (.go = 29.1473 nm"!) 
V'p = 3.250710-4 

'(40) 

The singular behavior of ^* (Ar*) and Xt (A''"*) of xenon 
can then be estimated by using Eqs. |[T]) to ((SJ, l|15p . and 
(|23|) , with xenon parameters of Eq. (|39|) . Such an esti- 
mation of dimensionless susceptibility (or dimensionless 
isothermal compressibility) will be represented by a full 
black curve (with label MR) in the next figures. 

As a most important result already underlined, d 
enables estimation of the extension of the Ising-Hke 
preasymptotic domain [see Eq. l(35|) ] 



•^PAD — 10 



(41) 



and the reduced crossover temperature [see Eq. ((36| ]: 



At 



A 



10" 



(42) 
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The temperature extension of the Ising-like preasymp- 
totic domain of xenon corresponds to the temperature 
range T—Tc < 30 mK [from Eq. l(4T|) ]. while the Ising-like 
predominant nature for crossover estimated by the mean 
crossover functions cannot extend beyond ~ Tc + 30 K 
[from Eq. (|42|) ]. More generally, the knowledge of -d also 
enables useful estimation of any relative temperature dis- 
tance and any effective value attached to a specific condi- 
tion of the mean crossover function on the complete At* 
range, as reported for example in lines 7 to 9 of Table [TTl 

Comparable values to the ones given by Eq. (|39l) can 
be found in several published papers but without explicit 
reference to the effective temperature range of the Ising- 
like preasymptotic domain. 

For example, in upper part (a) of Table IIIII , we 
have reported F"*" (column 3) and a+ (column 5) val- 
ues obtained by using different theoretical crossover 
functions calculated by several models labeled X = 
{GCrg4, MR6niax, MR, MSR,CPKFLM} (column 1) of 
respective Refs. (author?) [13, Efl, [13, 111 (column 9). 
These values result from fits of the isothermal suscepti- 
bility data of xen on p ublished by Giittinger and Cannell 
(GC) (author?) [2l|, with 7 (column 2) and A (column 
4) values fixed to their theoretical estimation. We note 
that the variations of the F"*" values are on a few per- 
cent level when the change in 7 values affects the third 
digit, while the variations of the a+ values cover a sig- 
nificative range such as 0.9 < 1.3, i.e., a+ ~ 1.1 
with ~ 20% deviation. Then, we have also reported in 
Table llllj the calculated values of the asymptotic scale 
factors (column 6), L^^^^ (column 7), and tpp (column 
8), by using Eqs. (|32l) to (|34l) . respectively, and an avail- 
able estimation of the leading amplitude = acS,^ of 
the correlation length (see column 7). We can observe 
a typical uncertainty of 25% in the 1} and ipp values. In 
addition, the theoretical error-bars and differences in the 
estimations of exponents 7, i^, A, and universal constants 
appearing in Eqs. ((32l) to (|34l) . can have comparable ef- 
fect on the estimation of t?, ipp, and go than the exper- 
imental uncertainties on the estimation of a+, F+, and 
^+ and the critical coordinates (such as Tc, pc, Pc, etc.). 

Nevertheless, accounting for the above variations of the 
1? values has no significant effect on the order of mag- 
nitude of ^pXb — 10"** estimated from Eqs. l(35|) and 
((39l) . We can then note that the apparent amplitude 
agreement arises in spite of the non overlap between the 
estimated temperature range of the Ising-like preasymp- 
totic domain [see Eq. I|4ip ] and the temperature range 
lO-'' < At* < 10-1 [i.e., 30mK < T - < 30 K], cov- 
ered by the Giittinger and Cannell's measurements (see 
also Appendix A). Therefore, the asymptotic amplitude 
evaluation provided by these fitting results cannot be eas- 
ily transformed in terms of the characteristic scale factors 
for two main reasons: 

i) measurements of the singular properties are made 
in a temperature range which do not reach the Ising-like 
preasymptotic domain; 

ii) fitting of the data is made with a Wegner-like expan- 
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Figure 1: Residuals -R% (kt)qq^ (expressed in %) between 
the dimensionless susceptibility data given in Table I of Ref. 
(author?) [2l| and the fitting expression of Eq. l|42|l : a) pa- 
rameters given in line GC4 of Table IIVI b) after the Tc shift 
of 0.5 mK to retrieve the 0.2% numerical level illustrated in 
Fig. 3 of Ref. (author?) [2ll|; In (a) are also illustrated the 
respectives deviations and related error-bars of the three "cal- 
ibrated" values of the isothermal compressibility for the three 
temperature distances indicated by the respective vertical ar- 
rows labeled At*^ , Atq2 and Ar^a (see text and Appendix). 



sion, whose validity is questionable outside the Ising-like 
preasymptotic domain (see I). 



Therefore, to replace with a suitable precision the in- 
dependent amplitudes of Eqs. ((39|) by the independent 
scale factors of Eqs. l(40|) . one needs to give the "rules" to 
interpolate from fitting results obtained in an experimen- 
tal range outside the Ising-like preasymptotic domain, to 
the hypothesized ones, only valid inside the Ising-like 
preasymptotic domain. 
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X 


7 


A 


r+ 


a+ 


Eq. I[32|l 


Eq. l|33ll 


Eq. dm) 


Ref 


GCrG4 

(pink) 


1.241 


0.496 


0.0577 


1.29 


0.0230 


24.305 
C(| = 0.184 nm 




(author?) [21| 


MR6n.a>. 


1.24194 


0.491 


0.057 


1.1844 


0.0194 


24.305 
= 0.184 nm 


3.25 X 10"'' 


(author?) [22| 


MR 


1.2395935 


0.50189 


0.057821 


1.23397 


0.02107 


25.6936 
= 0.184333 nm 


3.25 X 10"" 


this work 


MSR 


1.2396 


0.504 


0.0587 


1.11 


0.0171 


C(| = 0.184 nm 




(author?) [28] 


CPM 


1.239 


0.51 


n.a. 

(0.058)* 


1.08 


0.0162 


30.382 
= 0.184 nm (calculated) 




(author?) [26| 


LM 


1.240 


0.508 


0.0594 


n.a. 

(0.9)* 




n.a. 




(author?) [27] 



Table III: Amplitude-exponent results (colums 3 to 6) of the fit of the isothermal compressibility data obtained from the turbidity 
measurements of Giittinger and Cannell (author?) [2l|], along the critical isochore, using different theoretical crossover models 
(labels X, column 1) proposed in references given in the last column. Corresponding scale factor values (columns 7 to 9) 
calculated from Eqs. ll32ll to Ipj]! (see text for detail); n.a.: nonavailable. ; asterisk indicates a value used in the present work. 



# (fit) 


7 


A 


r+ 


a+ = a+ 






Tc 

(16.64°C) 


Ref 


GC/2 pink 


1.240 ± 0.002 


(0.496) 


0.0584 ± 0.0009 


1.07 ±0.06 


(0) 


(0) 


±0.5 mK 


(author?) \21] 


GC/4 pink 


1.246 ± 0.002 


(0.496) 


0.0551 ±0.0012 


1.62 ± 0.14 


-2.7 ±0.5 


3.6 ±0.8 


±0.5 mK 


(author?) \21] 


GC4 pink 


(1.241) 


(0.496) 


0.0577 ± 0.0001 


1.29 ±0.03 


-1.55 ±0.2 


1.9 ±0.5 


±0.5 mK 


(author?) \21] 


G3 


(1.240) 


(0.5) 


0.0574 


1.55 


-2.0 


(0) 




(author?) [12,1 



Table IV: Lines labeled GC/2, GC/4, GC4: Amplitude-exponent values obtained by Giittinger and Cannell from three rep- 
resentative fitting by Eq. ifiSjl of the isothermal compressibility data for xenon along the critical isochore; Fixed values of 
the parameters are given between brackets; Line labeled G3: Amplitude-exponent values obtained by a three-point calibration 
method proposed in Ref. (author?) [l3| and discussed in Appendix. 



3.2. Present status of crossover modeling in critical 

xenon 

3.2.1. Giittinger and Cannell 's analysis 



Due to the primary importance of Eq. l(32|) in a 
scheme where a single confluent ampUtude is readily ex- 
pected independent, in the Table IIV[ we have recalled 
the main informations given from Giittinger and Can- 
nell's experiment, which is the common support to com- 
pare theoretical flttings. Giittinger and Cannell observed 
that their susceptibility data systematically and contin- 
uously deviate from simple power-law behavior through- 
out the temperature range 9.6 x 10~^ < At* < 10"^ 
with the eflFective exponent increasing from 7e = 1.14 
to 7 = 1.246 ± 0.01. Thus the data appeared asymp- 
totically consistent with the theoretical estimations pro- 
vided at the end of the seventies, either using series 
calculation 7 = 1.250 ± 0.03 for the Ising-model (au- 
thor?) (stI. fHsl. or the renormalization-group result 7 = 
1.241 ± 0.002 for the Landau-Ginzburg- Wilson Hamilto- 
nian (author?) (HI, H^l- To support this conclusion, 
their susceptibility data were fltted to the following (four 



term) Wegner-like expansion 



K 



T,exp 



(At* 



r+ (At* 



^2x 



(At*) 



7 

2A 



1 



^3X 



(At*) 



3A 



(43) 



with A fixed at A = 0.496 (author?) [M] which was 
the theoretical value calculated at that time by the 
renormalization-group approach. In initial fittings, the 
parameters r+, 7, af^, aj^, a'^^, and Tc were indepen- 
dently adjusted, by variation of the data range in temper- 
ature. It was then shown that when the data range was 
narrowed to 9.6 x 10"^ < At* < 8.8 x 10~^ and the pa- 
rameter and a^^ removed (i.e., with a^^ — a'^^ — 0), 
the fit resulted in 7 = 1.240 ± 0.002 when 7 was freely 
adjusted, and the results for r+, af^ and Tc are in good 
agreement with all the other fits. This fit result is la- 
beled GC/2 and given in the corresponding line GC/2 of 
Table [IV] where the error bars quoted are one standard 
deviation allowing for the correlation between parame- 
ters, not the experimental uncertainty which increases 
for the data nearest T^, as discussed below. In fitting 
over the entire range, all the terms were required and the 
fit results (labeled GC/4) are given in line GC/4 of Ta- 
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ble llVI The Giittinger and Cannell's analysis to examine 
the consistency with the predictions by renormalization- 
group calculations was then made with 7 fixed at 1.241. 
The parameters for the fit result (labeled GC4) for the 
entire range in At* are reported in line GC4 of Table [iVl 



As noted by the authors, in all the fitting cases, the 
value of Tc found by the fitting procedure agreed to 
within 0.5 mK with the value observed by noting the 
temperature at which the meniscus formed upon cool- 
ing in small steps. It should be also noted that a shift 
of 0.5 mK in Tc amounts to a 2.5% change in fitted 
value at At* = 9.6 x 10~^, which was a substantial ef- 
fect in comparison to 0.2% deviations reported for all 
the fits (see Figure 3 in Ref. (author?) More 
generally, in the temperature range At* < 9.1 x 10~^, 
the authors have accounted for the large uncertainty in 
the corrections due to all the possible effects which are 
not at the level of about 0.2%, by using a shift in Tc 
in the fitting program as a mean of compensating the 
systematic errors which increase for the points nearest 
Tc. However, for the present discusion, it is essential to 
account for the real level of the fit deviation attached 
to the systematic experimental uncertainties. There- 
fore, we have illustrated in Figure [U the importance of 
the real uncertainty using the residuals (expressed in %) 

R% (K%-)nn ~ 100 ( 1 ) betwenn the raw data 

of the dimensionless susceptibility given in Table 1 of 
Ref. (author?) [2l[ and the calculated ones from Eq. 
(|43l) with the parameters given in Hne GC4. The part a 
of Fig. [H illustrates the 2% — 3% level of the true exper- 
imental error at At* = 9.6 x 10-^ while the part b of 
Fig. [1] illustrates similar residuals, but after the Tc shift 
of 0.5 mK, then reproducing Figure 3 of Ref. (author?) 
[m to show the 0.2% numerical level of the fit deviation. 



10= 10* lOf lOf lOi 



0.1^ I I I — I III 1 1 I I I — I |iii I I I I — I |iii I I I I — I |i 1 1 1 I I I — I pr 

o -7 /;:;.iooi 




10« 10' 10" 10° 10' 10^ 10^ 10' 1 



At 

Figure 2: (Color on line) Residuals R% (kt)gc^2 (expressed 
in %), as a function of Ar*, for the extrapolated descrip- 
tion of the isothermal compressibility of xenon within the 
Ising-like preasymptotic domain Ar* < £pad- Each curve 
label X refers to the line X in Table [TTH subscript indicates 
the number of terms of Wegner-like expansion (see also the 
text); Arrows indicate the different extensions of the Ising like 
preasymptotic domain (see text); horizontal (pink) segment: 
Giittinger and Cannell experimental range (author?) [2l||; 
open circle: temperature distance At*^ (see Fig. [TJ of the 
data calibration value; vertical (pink) segment: error bar of 
a calibrated valueof the isothermal compressibility of xenon 
from Hocken and Moldover measurements (author?) (HsI.IH^ 
(see text, Appendix A and Ref. (author?) [l^); upper hor- 
izontal axis: theoretical dimensionless correlation length cal- 
culated from the mean crossover function of Eq. IT]); line 
indexation with labels 1 to 8: je- or ^h-conditions reported 
in Table [n] (see text) 



In summary, comparing the three fit results of Table 
IIVI Giittinger and Cannell have shown that the xenon 
susceptibility was correctly represented by a Wegner-Hke 
expansion whose the two first terms may be interpreted 
as proposed by the renormalization-group theory with 
7 = 1.241 and A = 0.496 (author?) [60]. Despite a 
significative increase of experimental uncertainty as At* 
decreases, the uncertainty on the "mean" determination 
of the free leading amplitude r+ value can be estimated of 
the order of 2% — 3%, while the first confiuent amplitude 
aj value was found in the range 0.6 ^ aj < 1.45 and 
seriously affected by the presence of the other two con- 
fluent parameters and the fltted temperature range. The 
comparison between r+ and 0+ values given in Tables ITTTl 
and IIV| especially the low dispersion of the r+-values in 
a range 0.0570 < r+ < 0.0594, confirms an apparent re- 
duction of the a+-range, increasing the number of terms 
of the Wegner like expansion. 



3.2.2. Crossover analyses 

Since all the amplitude values given in Table IIIII origi- 
nate from the same experimental data, we consider that 
the extrapolation at At* ^ of the two-term expan- 
sion K^QCf2 obtained with the free parameters given in 
line GC/2 of Table HlH can act as an "ideal" experimen- 
tal result obtained within the Ising-like asymptotic do- 
main. Making then reference to this extrapolated exper- 
imental behavior, the relative comparison with the pre- 
dicted one using either a pure power law, or a two-term 
Wegner-Hke expansion, or a more extended Wegner-Hke 
expansion, or finally a complete crossover function, as- 
sociated to a: = {MR6max,MR,MSR, CPM,LM} mod- 
els (see Table IIIip . provides the required information 
to evaluate the Ising-like characterization in the ranges 
At* < £^Xd - 10"'' and At* > £^Xd - 10"^ sep- 
arately. We use the subscript 1, 2, etc., to refer for 
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one-, two-, etc., terms in the Wegner like expansion, 
while the abscence of this subscript indicates a complete 
crossover function (except for the MR6max case where 
the crossover function does not account for the classi- 
cal behavior close to the Gaussian fixed point). So that, 
anticipating the discussion given in Appendix, we have 
systematically enlighted (using a grey-blue area in the 
next figures), the temperature range in between the two 
temperatures T^Tc + lK (i.e., Ar*2 = 3.45137 x lO'^) 
and T = Tc-h 13.41 K (i.e., Ar^g = 4.62829 x lO'^) where 
the pVT measurements can provide two "standard" val- 
ues of the isothermal compressibility (see § A.3.b). These 
two temperatures are in the experimental range covered 
by the optical measurements of Giittinger and Cannell 
for At* > C^l^ ~ 10~*. 

Similarly, for At* < >CpXd — 10^^, we have also in- 
dicated in the following figures, the temperature T — 
Tc -t- ImK (i.e., Ar^j = 3.45137 x 10"^) "closest" to Tc, 
where a carefull analysis of the uncertainty on the theo- 
retical estimations of the Ising-like exponent 7 provides 
a cahbrated value of the isothermal compressibility (see 
below and § A. 3. a). This value is compatible with the in- 
terferometry measurements of Hocken and Moldover cov- 
ering the temperature range —bmK < T — Tc < 15 mK 
(i.e., -1.5 X 10-5 < At* < 5 X 10'^) (author?) 
As already shown, such a three point calibration can 
be used to define a practical three-term Wegner expan- 
sion with arbitrarily "fixed" values 7 = 77 = 1.240 and 
A = A/ = 0.5 for the exponents, which accounts for the 
isothermal compressibility of xenon in the experimental 
temperature range At* < AT^g, with ±1.5% uncertainty 
on the "online" values produced by the distinct measure- 
ment methods around these three temperature of cali- 
bration. The corresponding amplitude parameters of Eq. 
(|43|) are given in line G3 of Table iTVt and each associated 
fitting curve will be labeled G3 in the next figures. 

In Figure [21 we have reported the residuals (expressed 

in %), R%{k*j,)^^^^ ^ 100 ^ ^r^'^ -1^, as a func- 
tion of At*. Note that the indexation of the (lower 
and upper) horizontal axes of this figure are conform 
with the discussion of Table [ill while the experimental 
range of the GC measurements is given by the segment 
labeled GC where the open point indicates the temper- 
ature distance (noted At*j. in the following) where ligth 
scattering data are calibrated by the authors. The spe- 
cific analysis of the residuals obtained for the two-term 
expansion fitting (see all the curves labeled X2 in Fig. 
[21) demonstrates that the deviations in the temperature 
range £^Xd < ^'^* - 10"^- 10"^ induce a ± (1% - 2%) 
uncertainty level in the estimation of the F"*" value, which 
means that the value of the confiuent amplitude a+ is 
"undetermined" due to the correlation between free lead- 
ing amplitude and fixed theoretical value of 7. 

More generally. Figure [2] is a well-defined tool to show 
the relative importance of each term of the Wegner like 
expansion as a function of each decade variation of the 
theoretical correlation length (see upper axis). For ex- 



ample, the successive deviations at a 1% — 2% level of 
the curves Xi, X2, X, indicates that: 

i) the pure power law can be used in the range ^th ^ 

(1-2) X lo^ 

ii) the range £th > 1000 - 2000 needs to add the first 
confiuent term and effectively corresponds to the Ising- 
like preasymptotic domain, 

iii) a three-term-expansion seems sufficient to cover the 
range £th > 100 - 200, and, finally, 

iv) the range £th ^10 — 20 needs to use at least four 
or more terms in the Wegner-like expansion. 

To readily evaluate the relative importance of the ef- 
fects due to the 7 <^ r+ correlation, we note that the 
magnitude of the difference in 7 values here considered 
affects only the third digit of the theoretical estimations 
and cannot be related to the realistic error-bar (one or- 
der of magnitude greater) provided by the experimental 
determination of this asymptotic exponent (see also be- 
low the related insert of Fig. Therefore, in spite of 
the fact that measurements were performed beyond the 
Ising-like preasymptotic domain, the closest extrapolated 
Giittinger and Cannell's fitting to Tc remains highly cor- 
related to the theoretical value of 7, especially when Tc is 
fixed. This exponent-amplitude correlation was recently 
used by Luijten and Meyer (author?) [23| to re-evaluate 
the estimation of the leading amplitude for better agree- 
ment with the data for At* < 0.2 (see curve LMi in Fig. 
[21). However, such amplitude adjustment accounting for 
the data in a narrowed temperature range appears equiv- 
alent to the reverse effect of the theoretical min and max 
error-bar propagation outside the Ising-like preasymp- 
totic domain clearly shown in reference I (or in the Fig. 

of Appendix). Such an amplitude re-evaluation is 
then without gain on the Ising-hke asymptotic quality 
of the crossover analysis. For example, only by addition 
of the first term for confiuent corrections to scaling with 
a+ — 0.86 (see Table ITTTl fine LM), we can fit the exper- 
imental results (with deviations at a 0.5% level) on the 
complete experimental range (see curve LM2 in Fig. [21). 
The comparison of the curves LM2, GC/4, and GC4 in 
Fig. [21then shows that this "LM" fitting result which uses 
the parameters given in line LM of Table ITVl agrees with 
the initial four-term fitting results given in lines GC/4, 
and GC4 of Table [Tvl The respective ranges of the ampli- 
tude variations, r+ = 0.0584^°;°°;^° and a+ = 1.07j:|^:^^, 
illustrate the real difficulty to account for a single specific 
value of 7 in the range 7 = 1.24lg'QQ^. More generally, as 
illustrated previously in Fig. [H an "experimental" man- 
ner to account for each fixed value of 7 is to use Tc as a 
free parameter in fitting the GC's data. But, as observed 
by Giittinger and Cannell from their different fitting re- 
sults, the true uncertainty on the free parameter r+ is 
then at a few percent level. Finally, we can conclude that, 
due to the 2% uncertainty of the GC's measurements at 
the border of the Ising-like preasymptotic domain, the 
1% — 2% magnitude of the residuals for extrapolated the- 
oretical fittings at At* (< /^pXd) ~^ 0, is entirely related 
to the differences — 1.240 between respective leading 
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exponents, while the apparent reduction of the uncer- 
tainty in the aj-value results from the mandatory pres- 
ence of (at least two) supplementary terms to correctly 
account {~ 1%-level) for the experimental results in the 
temperature range At* > 10'^ > ^pXd — 10~^. 



3.2.3. Analysis with 7 fixed 

For a fixed theoretical value of 7 with Tc fixed, the 
leading amplitude r+ can be obtained only using addi- 
tional measurements performed well inside the Ising-like 
preasymptotic domain as defined by Eq. |[9]). The appli- 
cability of Eq. ([Ml) [equivalently, Eq. for the cor- 
relation length data] is then strictly limited within the 
Ising-like preasymptotic domain, thus requiring a pre- 
cision < 0.1%. In the absence of such an "ideal" ex- 
perimental result, one alternative way to eliminate the 
"amplitude-exponent" correlation on the leading power 
law term, is to introduce a smallest temperature distance 
to the critical temperature where a single experimental 
data for the isothermal compressibility acts as a standard 
value (author?) We note Ar^^ as this temperature 
distance such that Ar^j <C -CpAD (^^^ column 9, Table 
lll|) . For example, this method of calibration is detailed 
for xenon case in Appendix A, using the interferometric 
measurements of Hocken and Moldover (HM) (author?) 
[55.] in the range 2-3x10-6 < At* < 3-5 xlO^^ In Fig. 
[21 this experimental range is illustrated by the segment 
labeled HM. The selected smallest temperature distance 
is T - Tc = 1 mK, or Ar^^ (Xe) = 3.45 x 10"^, which ef- 
fectively appears within the Ising-like preasymptotic do- 
main, as evidenced by the corresponding vertical dashed 
line in Figure [21 Unfortunatly, the required precision is 
not reached and here we only discuss the relative com- 
parison between both (HM and GC) data measurements 
because there is no region of overlap between them and 
the two experiments rely on completely different effects. 
The interferometry measurements were performed in the 
very limited region of density and temperature quite close 
to the critical point such that the susceptibility data was 
represented by a pure power law with a highly correlated 
amplitude-exponent set. For example, the curve HMi in 
Figure [21 illustrates the deviation obtained using the re- 
sult 7 — 1.23 and r+ — 0.062 initialy found by Hocken 
and Moldover. The Hocken and Moldover measurements, 
reanalyzed by Sengers and Moldover (SM) with 7 = 1.24, 
yield r+ = 0.058, leading to the curve SMi in Figure 
[1 At Ar^j (Xe) = 3.45 x lO'^ the deviation of ~ 6% 
between HMi and SMi can be considered as representa- 
tive of the experimental uncertainty. These results seems 
to systematically deviate from fitting the raw data (full 
pink points) of Giittinger and Cannell's measurements. 
Moreover, accounting for the fitting results illustrated by 
the curves LM2 and/or GC/4, such a deviation appears 
to increase as the 7 value increases (or as the r+ value 
increases, equivalently). 

However, the Giittinger and Cannell's data are rela- 



tive to the isothermal compresibility value kt (Ar*^.) = 
1.9641 X 10-5 Pa"^ at Ar*, = 2.304 x IQ-^ (i.e., T^r = 
Tc + 0.6677 K) (author?) and the typical experi- 

mental uncertainty on the above reference value also is 
at a percent level. As a typical example, the difference 
on the GC value Pc (GC) — 1110 kgm^^ and our present 
value Pc = 1113. kgm"'^ (see Eq. l|A2p in Appendix) of 
the xenon critical density, contributes for a 0.6% differ- 
ence on the calibration. It is then interesting to evalu- 
ate the cumulative effect of a 0.5 mK shift of T^ and a 
1% change (uncertainty given by Giittinger and Cannell) 
of the reference value. The result is represented by the 
open pink points in Fig. [21 The fit deviation of these 
"corrected" points using our theoretical function is lower 
than 0.5%, as illustrated by the curve labeled MR in Fig. 

m 

Subsidiarily, Fig. [21 also shows that the available fit- 
ting results of the Giittinger and Cannell's data do not 
lead to the same standard value used to calibrate the 
mean crossover function. However, the attached relative 
difference appears compatible with an estimated error- 
bar of ±3% for this standard value, which is comparable 
to the fitting uncertainty attached to the leading ampli- 
tude due to the theoretical uncertainty on the 7 estima- 
tion. For example, considering the nine published pairs 
7 — (see Tables IIVI and the interferometry result 
of Refs (author?) [56.]) where 1.239 < 7 < 1.24194 and 
0.056 < r+ < 0.0594, the calculated mean value of the di- 
mensionless isothermal compressibility at T = Tc -I- 1 mK 
is «T,™ea„ = 3.4330 X 10^ with 3.3194 X 105 < k*^^^^^^ < 
3.5210 X 10^. Anticipating then the introduction of the 
calibrated value — 3.415 x 10^ discussed in Ap- 

pendix A, we note that the above mean value is 0.53% 
greater than the calibrated one, while the max and min 
values are 3.10% greater and 2.80% lower, respectively. 

As an undeniable consequence of the differences — 
1.240, (or the difference in leading amplitude, or the ab- 
sence of one precise data inside the Ising-Hke preasymp- 
totic domain, equivalently), the contribution of the 
highly correlated leading and first-confiuent terms propa- 
gates in any fit procedure which uses a complete Wegner- 
like expansion outside the Ising-like preasymptotic do- 
main, until a temperature distance of the order of At* ~ 
10"'^. We can then reasonably estimate (and precisely 
demonstrate below, in § 4.2 and Fig. [J]) that the fitting 
optimization occurs at finite temperature distance to the 
critical point (at least of the order of Ar* ~ 10~^), signif- 
icantly outside the Ising-like preasymptotic domain (and 
such as At* > At*^). In such a situation already known 
as due to the low convergence of the Wegner-Hke expan- 
sion, it is not easy to define with the needed precision a 
single crossover parameter which must account simulta- 
neously for two opposite roles: 

1) it controls the asymptotic universal features of the 
Ising-like singular behavior which has the restricted form 
of Eq. ([5]), but without constraint due to the lack of 
experimental data in the corresponding validity range; 

2) it detects the finite distance where classical-to- 
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critical behavior and/or non-critical behavior can occur, 
but may conjointly invalidate the use of the Wegner-like 
expansion with unique lowest confluent exponent. 

Therefore, in the absence of the "ideal" experiment 
within the Ising-like preasymptotic domain, the universal 
form of Eq. (|32l) can never be directly tested. This has 
three correlative consequences: 

i) the numerical value of d is undoubtedly related to 
the effects of many correction terms to scaling (for a more 
detailed illustration, see the discussion in terms of the ef- 
fective exponent and Fig. 4 of Ref . (author?) Q). To 
account for this essential result in the following, we note 

as the crossover parameter which is determined be- 
yond the Ising-like preasymptotic domain, to distinguish 
it from 1} of Eq. (fT3| which is the scale factor needed to 
characterize the asymptotic critical crossover within the 
Ising-like preasymptotic domain; 

ii) the use of Eqs. l(32|) to l|34p needs to have sup- 
plementary information such as, for example, the unam- 
biguous demonstration of the Ising-like uniqueness of the 
i^^-value (i.e. a constant z?£-value independent of the 
fitting domain and the fitted property); 

iii) alternative equations valid beyond the Ising-like 
preasymptotic domain which provide the equivalent scal- 
ing information contained in Eqs. l(32|) to (|34l) are suit- 
able to probe the identity d = dc and the application of 
the Hnear Eq. ^ for Ar* 0. 

We recall that the first attempt (author?) [121 to un- 
derstand the classical-to-critical crossover in xenon using 
the MR6niax scheme, was made from such a global anal- 
ysis of three singular properties (correlation length, sus- 
ceptibility and heat capacity) of xenon in a rather large fi- 
nite temperature range {T—Tc < 30 K, i.e. At* < 10^^). 
A rather well-defined value of = 0.0191 ± 0.0095 was 
then obtained, which was subsequently analyzed in terms 
of the infiuence of the corrections to scaling in real sys- 
tems, even outside the Ising-like preasymptotic domain. 
The small difference with the present hypothesized value 
d = 0.021069 [see Eq. gO])], only refiects the error-bar 
propagation of theoretical exponent and amplitude esti- 
mations in the respective fitting of the MRGmax or MR 
discretized values, and not the true experimental uncer- 
tainty (see a detailed discussion of the experimental un- 
certainty in Appendix A). In other words, the hypothe- 
sized asymptotic value of d appears independent of the 
Ising-like Eq. (|32l) . Similarly, in the other crossover mod- 
eling, the uniqueness of the crossover temperature scale 
appears independent of the theoretical approach in the 
sense where it is de facto obtained at a "single" finite 
temperature distance, noted (Ar*)^^. Here single means 
that, in a log-log scale, the value of the free parameter 
is optimized on a fine temperature domain (covering typ- 
ically less than one decade) centered on a large temper- 
ature distance (At*)^^, such that ^ ^^♦''^ ^ 1- For ex- 
ample considering the Giittinger and Cannell's measure- 
ments in the next Section, we show clearly that (At*)^^ 
corresponds approximately to the "largest" experimental 



value (At* 



0.05 - 0.1 < At^ of the fitted temper- 



ature range reported in Fig. [2l Then the ratio ^ a^*''^ 

is greater than 10''. This value is definitively outside the 
Ising-like preasymptotic domain and similar to the value 



of crossover temperature scale At^ 



M 



0.23 obtained 



from fitting using CPM and MSR crossover descriptions. 

The above analysis Hmited to the results of Figure 
[2] extrapolated within the Ising-like preasymptotic do- 
main, justify our special attention in the next section to 
methods which account for the calibrated contribution of 
the theoretical leading power law, in order to check in a 
self-consistent manner the iJ^-determination beyond the 
Ising-like preasymptotic domain. 



4. XENON CHARACTERIZATION BEYOND 
THE ISING-LIKE PREASYMPTOTIC DOMAIN 



4.1. Ising-like nature of the dimensionless scale 
factors 



The finite and restricted At* range where the mean 
crossover functions correctly represent and ^exp 

data, determines the effective extension £ead the 
Ising-like extended asymptotic domain (EAD). Since the 
theoretical expressions of Eqs. ^ and ([2]), are in 
the form of a complete crossover, -Cead involve 
correction-to-scaling terms higher than the first one. In 
such a situation, the following condition >CpAD < ^ 
'^EAD occurs. This is precisely the case for the isothermal 
compressibility data of xenon obtained by Giittinger and 
Cannell. As a consequence, the value of i?£ introduced 
by the analytic relation t — i^^At*, is related to an "un- 
defined" domain of extension -Cead- Here, "undefined" 
means "beyond the Ising-like preasymptotic domain", so 
that we are not able to appreciate the effective infiuence 
of all the numerous corrections neglected in the massive 
renormalization scheme of the <i>^ model (see I) . We must 
solve new correlative difficulties concerning the effective 
number (which can thus be greater than 3) and the na- 
ture (which can originate from the neglected analytical 
and confluent effects in the critical massive renormaliza- 
tion scheme) of the fluid-dependent parameters. 

Indeed, to complete Tc and ac knowledge in the ab- 
sence of information concerning the true range of the 
Ising-like preasymptotic domain for an actual fluid, it 
is proposed (author?) 0, [l3| to replace tp, and go 
by three new adjustable parameters Xg £, Lq £, and 'ffc, 
modifying Eqs. (|23|) and l(20|) in the following way: 



„(Ar-) 



(At*)^ 



(L*J-^Z+(At*)^ 



(44) 



(45) 
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with 

t = ??£Ar* 
The new characteristic set 



I (/3c) ^;ac;Sicpx} 



(46) 



(47) 



must substitute the set Qc (Ar* < ^ead) of Eq. (|25l) . 
while the new dimensionless set 



Si 



CP,C 



(48) 



must substitute the set Eisf of Eq. l(26|) . The subscript 
ICP, C recalls for a single crossover parameter to charac- 
terize the crossover behavior observed in the finite tem- 
perature range >CpAD < ^'^* ^ ''^ead- 

In comparison to Eqs. l(23|) and l(20|) . the noticeable 
modification of Eqs. (|44|) and (|45| is the leading term 
in which is no longer involved in the asymptotic scal- 
ing part of the critical behavior expressed in terms of 
the physical field Ar*. As introduced, Xq ^ and Lq ^ are 
prefactors for each corresponding property, here selected 
as independent and characteristic of the fluid by virtue of 
the two-scale-factor universality (provided that the same 
length unit was used to define the dimensionless quanti- 
ties (author?) [2^). Correlatively, i?£ is a, pure crossover 
parameter, with same value above and below T^, which 
exclusively controls the magnitude of many correction 
terms to scaling. In addition, 'ffc can also integrate some 
effects of the neglected terms Hnked to the supplemen- 
tary confluent exponents, such as A2 or A3, accounting 
for practical numerical approximations such as A2 « 2A 
or A3 « 3 A, or the effects of using T to replace Tc in the 
energy unit and in the dimensionless form of the tem- 
perature distance to the critical temperature. For a fluid 
/, the determination of i9c is then equivalent to the de- 
termination of 'Cg^Q. In this general case, the physical 
leading amplitudes can be calculated using: 



(49) 



(50) 



i.e., without reference to -de- However, the subscript C 
recalls for the determination of ■)?£, then correlatively, 
Lq £ and Xq £, made beyond the Ising-like preasymptotic 
domain. The single (system-dependent) first confiuent 
amplitude can be calculated from '&c uniquely, using one 
independent equation among the two following equations: 



(51) 



As suggested in Ref. (author?) 0, from similar fit- 
ting of the correlation length and the susceptibility in the 



inhomogeneous domain, the specific heat in the homoge- 
neous and non homogeneous domains, and the coexisting 
density measurements in the inhomogeneous domain, one 
must verify the uniqueness of the value (along the crit- 
ical isochore). Considering several properties allows then 
consistent tests for the determination of the set Sicp,c^ 
in coherence with the basic hypotheses of the renormal- 
ization group approach at the origin of the theoretical 
crossover functions. 

To avoid this large task, in a first approach, we con- 
sider xenon for which we hypothetize the existence of a 
single scale factor '&c in the temperature range Cp^^^ < 
At* < >Cead- Fitting then and ^*^p data in this 

temperature range with Eqs. I|15p and l|23p . produces the 
set Sicp,Ct i-e-, the determination oi'&c- Extrapolating 
the results for At* 0, we can identify the value of the 
asymptotic scale factor to 



13 = dc for Ar* < C^lj^ 



(53) 



After that identification, we can introduce the dimen- 
sional prefactor Lqx = acLg £ leading to define the 
asymptotic wave number go of Eq. (|33l) by the follow- 
ing parameter: 



50 = gox 



(54) 



The two remaining asymptotic scale factors L^^-''^ [see 
Eq. IpT]) ]. and Vp [see Eq. (flSl) ]. are then obtained by 
the hierarchical equations: 



acgox 



{1/} 



(55) 



(56) 



Finally, for Ar* 0, the scale factor set 
SsFX = {^c;^^'^^;^px} with < Ar* < (57) 

has the appropriate asymptotic form to compare with 
SsF [Eq. ((26l) ]. except the noticeable subscript £ which 
recalls the "non-asymptotic" Ising-Hke nature of these 
numbers which originates de facto from the determina- 
tion of dc, Lg £ and Xg £ made beyond the Ising-like 
preasymptotic domain. 

However, when this "hypothetical" xenon is character- 
ized by Eq. (07]) [or Eq. ([57])], thus Eq. ([55]) is true and 
the following variable transformations 



ac 
At* 



4h = L^i/>r 

t = ^cAt* 



(58) 
(59) 



result in the asymptotic collapse (over the extension 
(52) ^PAD (author?) [l3,[ll|) of any physical curves of equa- 



tion ^ (Ar*) into the universal curve of equation 

4h (t) = L^i/>r (Ar*) = goxUv (Ar*) (60) 
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We have used this universal feature due to the scale 
factor nature of i?£ to obtain the xenon values of lines 8 
and 9 of Table ITTl from theoretical values of lines 5 and 6, 
respectively. In such a numerical approach, the "hypoth- 
esized" xenon values of the characteristic set Sicp,c are 
the following [see Eq. IjAlop in Appendix] 

f = 0.021069 1 

Sicpx = I K.c = 0.443526 i (61) 
[ X*^ = 0.214492 j 

with 

Lox = acL5,£ = 0.390972 nm (62) 
and, obviously 

h^^f^ = ttcgox = 25.6936 (63) 

More generally, our hypothesized values of Eq. l(6T|) pro- 
vide the (expected) identity Ssf.,c = Ssf- Accordingly, 
the universal scaHng form of the correlation length of 
xenon reads as follows 

L{i/}r(Ar*) = ^^^gg^^4h(t) (64) 

[mn] 

with 



For example, looking back on our previous analysis 
of the residuals for the xenon isothermal compressibility 
case, by reversing Eq. l(64|) . we can transform At* in uni- 
versal values of £th, as illustrated in the upper axis of Fig- 
ure [21 Simultaneously, as previously mentioned, we are 
also able to illustrate in this figure (or in any figure which 
use a xenon At* -coordinate) all the conditions estimated 
in Table [III In the part of the Giittinger and Cannell's 
experimental range which corresponds to the theoreti- 
cal condition 4h ^ (2.5 — 3) L^^^^ ~ 70, the universal 
behavior of the functional form x* (4h) can then be an- 
alytically "tested" in this Ising-like extended asymptotic 
domain. Simultaneously, we will also confirm below the 
uniqueness of the crossover parameter by using an alter- 
native facet of the universality accounted for by the mean 
crossover functions. As a matter of fact, to illustrate the 
asymptotic Ising-like transformation of each thermody- 
namic property P* (At*) attached to the uniqueness of 
the crossover parameter, we can also use the effective 
universal behavior of the local exponent ep^th (t) to re- 
place the one of the correlation length £th (i) , in order to 
construct the "universal" scaHng form Fp (ep_th)- There- 
fore, in a second approach with the objective to illustrate 
this transformation in a self-consistent manner for the 
susceptibility case, we consider the universal and exper- 
imental effective amplitudes attached to the local power 
laws with effective exponents, such as introduced in Ref. 
(author?) J^. 



t = 0.021069AT* 



(65) 



Either the "theoretical" value 0.0343085 nm (= [gox] ^) 
of the length unit in Eq. l(64|) , or the "measured" value 
0.184333 nm (= S,qc) of the correlation length amplitude 
in Eq. 128]) , cannot easily be related to a real microscopic 
length of xenon atom, while the fact that the "fitting" 
value 0.390972 nm (= Lqx) of the dimensional prefac- 
tor is comparable to the size of the xenon atom can be 
considered as a fortuitous result (author?) [HH . 

As expected, more fundamental is the dimensionless 
nature of Eq. ([631 which can then provide the needed 
uniqueness of the length unit for better understanding 
of universality. Indeed, the above asymptotic collapse of 
the correlation length implies equivalent asymptotic col- 
lapse of any other singular thermodynamic property, by 
virtue of hyperscahng. That means that, when the scale 
factors L^^^^ and are determined in the Ising-like ex- 
tended asymptotic domain ^pXd < At* < >Cead' the 
validity range where the universal collapse onto the the- 
oretical thermodynamic behavior is expected, can be al- 
ternatively discussed in term of the value £th {t) > L''^^-'^-'^. 
In other words, all the crossover functions are "Ising-like" 
universal: 

i) only over the temperature range where the crossover 
parameter is unique; 

ii) only for a single dimensionless critical length L^^-'''^ 
common for all the fluids which obey to this single pa- 
rameter crossover description. 



4.2. Effective exponent and amplitude beyond the 
Ising-like PAD 

>From Xth (t) of Eq. the related local values of 
the effective exponent 7e.th {t) and effective amplitude 
Z+ g (t) are given by the equations 



7e,th (t) = - 



dLn [xth jt)] 
dLnt 



Xth jt) 



(66) 



(67) 



Eliminating t [then simultaneously ehminating the scale 
factor z?£, since t = 'i?£AT*], the theoretical classical-to- 
critical crossover is characterized by a "universal" curve 
^J,e (7e,th) over the complete range 7mf < 7e,th {t) < 7 
(see curve labeled $3 (l)-MR in Fig. [3|). 

Our present interest is restricted to the Ising-like range 
7e,th {t) > 71 = ^"'"2'"'' (see the corresponding Ising-like 
range defined in the upper part of Fig. [3|). The theoret- 
ical Ising-like limiting point takes universal coordinates 

I7; ^1 (upper cross in Fig. [3|). The "small" exten- 

sion 7-7e,th < ^x+A (/:p7]!)'^ ~ 0.006 of the Ising-like 
preasymptotic domain is magnified by the curve labeled 
$3 (l)-MR in the insert of Fig. [3 On the other hand. 
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^3(1) - MR 
Xe - MR 
Xe - GC4 
Xe-pVT 



1 1.05 1.10 1.15 1.20 1.25 
% 1.234 



Figure 3: (Color on line) Two-parameter transformations f^^ {^q,Cj'^c) (see arrays), at 7e = constant value, from the theoretical 
universal behavior of the susceptibility (upper curve), to the xenon effective behavior of the isothermal compressibility (lower 
curve) (in dimensionless variables); Xq £: adjustable parameter for the general case of Eq. I|44p . unequivocally defined from 
comparison in positions of the two Ising limiting points (right crosses); adjustable parameter unequivocally defined by the 
two other points, either within the Ising-like preasymptotic domain (see the insert and the text), or outside the Ising-like 
preasymptotic domain (full black circles T and P on the two curves labeled 'li's (1) — MR and Xe — MR, respectively). Note 
the inversion of the relative position between the curve labeled $3 (1) — MR and the curve labeled ar which is related to the 
asymptotic (logarithmic) approximation of Eq. I|68p . beyond the Ising-like preasymptotic domain. 



the curve ax corresponds to the asymptotic behavior of 



the derivative 



7c ^7 



of equation 



7-7e,th 

A|zi-+| 



1 - log 



7-7e,th 



(68) 

The vertical double arrow with label(l)indicates the 
above logarithmic divergence of (^r^j . We note 

the significant difference between the curve ax and the 
curve {S) which results from "analytic" error-bar corre- 
lation between the Ising values of 7 and (Zj) . As a 



matter of fact, the curve (5) corresponds to the linearized 



slope C+o 7„i„-7. 



between the respective 



bounded coordinates of points A and B (see inserted ta- 
ble in Fig. [Sand Ref. (author?) for data sources). 

Since only two parameters (Xg ^ and dc) are free in 
fitting Eq. l(44|) . Fig. [3] illustrates how the adjustable 
prefactor Xq ^ acts at the exact value of the Ising ex- 
ponent, since each Ising point of the fluid takes the co- 
ordinates F'^ — Xg^£ (^x) represented by the 
lower cross in Fig. [3] for the critical xenon case. The 
resulting prefactor (Xg = Xg £ is schematized by 
the double array between two crosses in Fig. [31 There- 
fore, caHbrating the r+ value as suggested in previous 
paragraph, fixes this prefactor value which governs the 
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universal collapse of the Ising-like limiting point through 
Eq. l(50l) . Xq £ acts in a equivalent manner to L^^^^ 
[see Eq. (|64|) ] for the correlation length case. In prin- 
ciple, the critical divergence in the initial slope at the 
limiting points provides the second "Ising-like constraint" 
which is needed to determine the asymptotic scale factor 
i). However, our previous description of the Ising-like 
preasymptotic domain and its above geometrical illus- 
tration underline the challenging (theoretical and exper- 
imental) difficulties to provide the exact characterization 
of the asymptotic scaling when a property reaches the 
Ising-like limiting point along a curve of "universal", but 
infinite, slope for all the physical systems. 

At contrario, the description of the 7e-variation in the 
range 7 — 7e ?J 0.015 (i.e., in a temperature range sig- 
nificantly beyond the Ising-like preasymptotic domain), 
appears now simplified, using precisely the Giittinger and 
Cannell's results of Eq. (|43l) to define the following ef- 
fective exponent by : 



7e,exp (At*) = -- 



dLr 



(At*)] 



dLn (At*) 



and its attached effective amplitude by : 



r+(Ar*) 



T,exp 



(At*) 



(At* 



(69) 



(70) 



The resulting single curve r+ (7e,exp) is illustrated in 
the lower part of Fig. [3] (see curve Xe — GC4), while the 
expected (two parameter) transformation /(Xq^jI?^) 
able to insure the universal collapse between the theoret- 
ical curves (labeled Xe — MR and $3 (l)-MR), is schema- 
tized by a double array between two points on each curve 
of well-defined finite slope. This transformation must 
contain the needed both constraints for the (point) po- 
sition and the related (tangent) direction. Therefore, 
the scaling nature of the collapse beyond the Ising-like 
preasymptotic domain is significantly different in fitting 
procedure which either eliminates, or accounts for, the 
contribution of the leading term. In the latter situation, 
we can then replace the prefactor Xg ^ by the calibrated 
leading amplitude F"*", as seen above. 

In the first case, at large temperature distance, the fit 
procedure based on Eq. (|44l) is mainly equivalent to a 
dominant constraint in "direction" given by the following 
relation between the two effective exponents: 

7e,exp(AT*) = 7e,th[^?£(AT*)] (71) 

We numerically solve Eq. l(7T|). using the Giittinger and 
Cannell's fitting results given by Eq. l(43|) . then providing 
the 7e [(At*)] and -de (At*) values as a function of At*. 
Both results are reported as the curve labeled GC4 in 
Fig. UK [7e as a function of At*], and the curve labeled 
1 in Fig. [Hd [dc as a function of At*], respectively. 

In the second case, to account for the contribution of 
the leading term needs to use the following scaling rela- 
tion between the two effective amplitudes 



(72) 



Now, the transformation / (Xg 'dc) = i'&c. 



,7-7e 



IS 



explicit in Eq. (|72l) . Its takes a convenient effective power 
law of the crossover parameter while the prefactor 
Xg £ has (as expected above) the same value whatever the 
7e ( = 7e,exp — 7e,th) value is. The constrained "position 
and direction" are accounted for correctly. Therefore, 
using Eq. ((50l) to eliminate Xq£, infers the pure dc- 
dependence of equation 



F+ 

r+ 



7-7e 



(73) 



By appropriate combination between the Giittinger and 
Cannell's fitting results and the mean crossover function 
for susceptibility, we numerically calculate the local value 
of the crossover parameter over the complete experimen- 
tal temperature range, using the equation 



1 



F+ 

^ e 

F+ 



(74) 



It is essential to note that for each fiuid for which r+ is 
known, Eq. I|74p . applied in the extended asymptotic do- 
main At* < >Cead ' takes equivalent "Ising-like" meaning 
to Eq. ((32I) applied within the Ising-like preasymptotic 
domain At* < -CpAD- In Figure Hb, we have reported as 
a curve labeled 2, the calculated value of 'dc as a func- 
tion of At*, using Eq. ^ with r+ = 0.0578204 [see Eq. 
([39l) ]. The available part of these curves 1 and 2 must 
be restricted to the experimental temperature range il- 
lustrated by the segment labeled GC. 

For both cases, the t^^-change as a function of At* 
within the range lO'^ < At* < lO'^ of Fig. [4b, can 
be approximated by our hypothetical Ising-like asymp- 
totic value 1? = 0.021069 of Eq. (091). In Fig. Hh the 
corresponding curves (labeled 1 and 2) of the residu- 
als R% {dc)^ = 100 - 1) (expressed in %), together 
with their mean curve (labeled m), are given. In the tem- 
perature range 5 x 10^^ < At* < 10^^ of Fig. H):, the 
green area corresponds to the error-bar of ±15% for both 
determinations. Even at such large values of At*, the 
main significant result obtained from the massive renor- 
malization scheme, is the estimation of the temperature- 
like crossover parameter in conformity with the asymp- 
totic two-scale-factor characterization of the fiuid, espe- 
cially using the true scaling Eq. l(74|) . 

Such a temperature range where a "measurable" value 
of the exponent difference 7— 7e occurs, was largely inves- 
tigated in the seventy's (author?) [12], when the scaling 
approach of the fluid universality was based on the ef- 
fective "universal" values of the critical exponents - as 
for example 7e = 1.211 (author?) p3], "feos = 1-19 
(author?) [Ij, Si], and 7e = 1.16665 (author?) [H] 
- involved in effective power laws and/or effective "uni- 
versal" form of a rescaled equation of state. Anticipating 
a more detailed discussion given in Appendix A, we can 
use the data reported on Table [V] for these 7e values 
to easily demonstrate, using the corresponding arrows in 
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10'= lO-'^ 10-^ 10-^ 10-^ 10-1 1 101 10=* 

At 

Figure 4: (Color on line) Ising-like asymptotic behaviors of: 
(a), the effective exponent 7e for the xenon isothermal sus- 
ceptibility, (b), the xenon crossover parameter -ffc, and (c), 
the percent deviation between i?£ and the xenon asymptotic 
scale factor calculated from the scale dilatation method, as 
a function of the relative temperature At* (lower horizontal 
axis) and the theoretical correlation length £th (upper hori- 
zontal axis); curve labeled MR: from MR crossover function 
of Eq. I|23p and xenon parameters of Eq. I|40p : curve GC4: 
from a fit of the Giittinger and Cannell's data by Eq. I|43p 
and xenon parameters in Table [W] line GC4; curve 1: from 
Eq. I|7ip : curve 2: from Eq. ([74]); Curve m: mean value of 
curves 1 and 2; lines labeled CPM and MSR: with i9x-values 
of Table Hill obtained from the fitting results of the Giittinger 
and Cannell's data by CPM and MSR crossover functions (see 
text). 



Fig. SK, the "Ising-like nature" of the covered temper- 
ature range At* < 0.05 — 0.1 ^ Ar^, or alternatively 
but equivalently, the "Ising-like nature" of the covered 
correlation length range £th ^ 70 > 2.5L^^-^h That also 
gives interest to revisit (author?) the effective uni- 
versal formulation of a parametric equation of state using 
the master crossover functions to validate the universal 
features observed in the well-defined Ising-like extended 
asymptotic domain of the fluid subclass. 

Finally we note that, when Tc and ac are known, the 



mean crossover functions take a convenient controlled 
form to determine a single crossover parameter -de in 
a temperature range beyond the Ising-like preasymp- 
totic domain. The Ising-like nature of this crossover pa- 
rameter is then revealed by the dimensionless value of 
a single characteristic length. However, to define the 
"minimal" set of three characteristic parameters such 
as {"(?£; Xq £; ^j, or, alternatively but equivalently, 
{d;ipp; go}, and {a+; r+; ^^J"}, one needs a "data calibra- 
tion" from measurements performed within the Ising-like 
preasymptotic domain. The practical interest of this data 
caHbration is given in the Appendix A, without reference 
to the estimated precision of the experimental method 
able to provide such measurements in this closest tem- 
perature range. 

In the absence of explicit thermodynamic definition for 
the prefactors Lq £ , £ and the crossover parameter 
■dc [or for the scale factors go, 1}, and ipp], the remain- 
ing difficulty is to compare between distinct fiuids which 
show differences in their fluid-dependent amplitudes , 
r+, and a+. This difficulty can be solved by applica- 
tion of the scale dilatation method when the localization 
of their liquid-vapor critical point is known, as shown in 
Appendix A for the xenon case. 



5. CONCLUSION 

Using xenon as a standard critical fiuid, and the mean 
crossover function for susceptibility in the homogeneous 
phase as an illustrative example, we have estimated the 
values of the fiuid-dependent parameters which are com- 
patible with the universal features predicted by the mas- 
sive renormalization scheme. A special mention for the 
three dimensionless parameter characterization within 
the Ising-like preasymptotic domain was given, in spite 
of the large theoretical and experimental uncertainties at 
such "closest" temperature distance to the critical tem- 
perature. Using the abundant literature now available 
from several fittings of Giittinger and Cannell's data, 
we have demonstrated the great advantage of the mean 
crossover functions to provide an unambiguous determi- 
nation of one fiuid-dependent crossover parameter out- 
side the Ising-like preasymptotic domain. Specifically, we 
have clearly shown that the value of this crossover param- 
eter is entirely governed by fitting the data at "largest" 
distance to the critical point, leading for example to an 
apparent reduction of the uncertainty in the determina- 
tion of the amplitude of the first confiuent correction- 
to-scaling term. Finally, the magnitude of the resulting 
deviations and the range of temperature where it is to be 
observed are exactly accounted for. 

However, in the absence of controlled information 
about the Ising-like preasymptotic domain description, 
only the similar use of the complete set of mean crossover 
functions applied to several properties is able to demon- 
strate that the crossover parameter conserves its "Ising- 
like" nature, even outside the Ising-like preasymptotic 
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domain. Alternatively, it was recently proposed in II 
and herafter justified in the Appendix A for the xenon 
case, an appropriate modification of the mean crossover 
functions which only uses three master (i.e. constant) 
parameters. The modified crossover functions represent 
the master (i.e. unique) singular behaviors of the one 
component fluid subclass in a well-deflned Ising-like ex- 
tended asymptotic domain. In such a situation, the real 
extension and amplitude of the singular behavior of the 
fluid properties can be estimated for any one-component 
fluid for which the vapor-Hquid critical point is locaHzed 
in the pVT phase surface. 
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Appendix A: ISOTHERMAL COMPRESSIBILITY 
OF XENON 

In the flrst part of this Appendix we calculate the char- 
acteristic parameters of xenon involved in Eqs. (|38| and 
(l39l), b y using the scale dilatation method (author?) 
[la . fl4l | . The needed information is then given by the set 
Q™™ made of four critical coordinates of the xenon crit- 
ical point. Indeed, for any one-component fluid, Q™" 
reads as follows (II) 

QTZ = {Tc;Pc;vpy,j',} (Ai) 

where vp^c = ^ is the critical molecular volume and 7^, is 
the common critical direction at the critical temperature 
of the critical isochoric line and the saturation pressure 
curve in the p; T diagram. 

The second part shows that the xenon parameters 
given by Eq. (|40l) can be used as entry data for the 
Eqs. ((3]) and l|23p to represent accurately the singular 
behavior of the xenon isothermal compressibility in the 
temperature range such as T — Tc < 5 — 10 K. 

The last part discusses the data calibration and the 
uniqueness of the Ising-like crossover parameter within 
the extended asymptotic domain. 



1. Xenon data sources from application of the scale 
dilatation method 

a. Xenon critical coordinates 

The selected critical coordinates of xenon are: 

Tc = 289.733 K; 

Pc = 5.84 MPa; 

Pc = 1113 kgm-3 or Wp^c = 0.19596 nm3; (A2) 

7^ = 0.1192 MPaK"\ 

withTOp = 2.181 X 10-25 kg. 



They result from the combined analysis (author?) 

of pVT measurements of Habgood and Schnei- 
der (author?) ji^, and coexisting density measurements 
of Weinberger and Schneider (author?) |f45|, Cornfeld 
and Carr (author?) {s^l, Thoen and Garland (author?) 
[67| . and Balzarini and coworkers (author?) |70l. l7ll|. 

The xenon critical temperature and pressure were flxed 
to the values recently provided by Gillis et al (author?) 
[ll which are compatible with = 289.740 ± 0.003 K 
and Pc = 5.8400 ± 0.0005 MPa obtained from Schneider 
et al's measurements. 

The Pc value has an uncertainty of ±5kgm-^ (~ 
0.5%), which accounts for the pc values of Schneider 
et al's {pc = 1110±2kgm~3) (author?) [H, 
Cornfeld and Carr's {pc = 1111.2+|;4kgm-^ for three 
different estimations) (author?) [63], Baidakov et al's 
{pc — 1112. 8±n. a. kgm"^) (author?) [64], and Balzarini 
et al's {pc = 1099 ± ri.a.kgm-3 (author?) [70], Pc = 
1116.0 ± 1.7kgm-3 and pc = 1114.0 ± 1.7kgm~3 (j^^. 
thor?) 

The value 7^ = 0.1192 ± 0.0005 MPa K"^ (~ 0.5%) 
was estimated by one of us (author?) |12] by graphi- 
cal analysis of the pVT measurements of Habgood and 
Schneider, which agrees to other literature values 7^, = 
0.1192 MPaK"^ (author?) Q, 7^ = 0.1196 MPaK"^ 
(author?) [H3, and 7^ = 0.120 MPaK"^ (author?) 
[eil. The -f'c value used in Ref. (author?) ^ differs 
by 0.19%. 

More generally, we note the remarkable agreement with 
the critical set defined by Gillis et al (author?) [t^I in 
their recent analysis of the sound attenuation (in the fre- 
quency range 100 < / (Hz) < 7500) within thermoacous- 
tic layers between solid surfaces and xenon at critical 
density. 



b. Physical and master amplitudes from the scale dilatation 
method 

>From Eq. (|A2|) . the critical values [see Eq. ([Ml)] of 
the energy (/3c) and length ac units of xenon are the 
following: 

(/3J-1 = keTc = 4.0003 x lO'^i J 

/k T\i (A3) 
ac=(!^j = 0.881508 nm 

while the values of two xenon scale factors Yc and Zc are 
the following: 

Zc = -^jP^ = 0.28601 , . , 

Yc = jc— - 1 = 4.91373 ^ ^ 

We have calculated the values [see Eq. l(39l) ] of the xenon 
amplitudes r+, , and aj, (or a^) by using the follow- 



20 



ing relations : 



r+ = 



{ycY 



= 1.23399 
= 0.209111 



= 0.0578204 



(A5) 



(with ac^+ = 0.184333 nm) 
ective v 

-E+, and (or Z^'^) are (author?) (Tj 



where the respective values of the master amplitudes Zt , 



xenon reported by Giittinger and Cannell (author?) 

we note that the turbidity data are fitted in the 
Orstein-Zernike approximation within a 1% precision in 
the temperature range 28 mK < T — < 2.54 K [i.e., 
9.6 X 10^5 < At* < 8.76 x lO^^j ^ ^gj^g the present the- 
oretical estimation of the correlation length and isother- 
mal compressibility. In particular, we estimate the ref- 
erence value T (Tor = + 0.6677 K) = 4.1067 m"! of the 
turbidity, in excellent agreement with the Giittinger and 
Cannell's one t {Tg^) = 4.1 m~^ (author?) [2l|, without 
any adjustable parameter. 



2+ - 
7+ — 



0.570481 

0.119 

0.555 



(A6) 



Hereabove, universal features within the Ising-like 
preasymptotic domain are correctly accounted for by the 
following equations 



4" 

with 



0.37695 



% = 0.67919 



0.83812 



(A7) 



The amplitude value £,q ~ 0.184333 nm estimated using 
the scale dilatation method, compares favorably with the 
one = 0.1866 ± 0.001 nm recently used by Gillis et al 
(author?) [t^I to analyze sound attenuation within ther- 
moacoustic layers between soHd surfaces and xenon at 
critical density. Moreover, this amplitude value is also in 
agreement with the following ones obtained from analyses 
of (static and dynamic) light scattering measurements: 
(i) £_Q — 0.2 nm, with v = 0.63, in the temperature range 
22 mK < T-Tc < 3.3 Kli.e,j^8x lO""^ < At* < IQ-^] (see 
Refs. (author?) 0, Hfl, il) ; (ii) = 0.1934 nm, with 
v = 0.62, in the temperature range 2.6 mK < T — Tc < 
10 K [i.e., 9 X 10-^ < At* < 3.4 x lO-^j from Ref. (au- 



thor?) [Sai; (iii) ^0 = 0.184 ± 0.009 nm and one conflu- 
ent correction term (with — 0.55 and A = 0.5), in 
the temperature range 28 mK < T — Tc < 3.65 K [i.e., 
9.6 X 10-5 < At* < 1.26 x lO'^] from Ref. (author?) 

Using the above xenon parameters in the theoretical 
estimation of the correlation length, we also underhne 
that the agreement with the experimental measurements 
extends to the range ^ ^ 3, i.e., in a temperature range 
At* < L^%jj ~ (2 - 3) X 10"^ which extends largely be- 
yond the Ising-like preasymptotic domain (see Ref. (au- 
thor?) [i3|). For example, at the calibration temper- 
ature T = Tc + IK (At52 = 3.45137 x lO'^, see be- 
low), our calculated value of the correlation length is 
^ — 68.6323 A (i.e., ~ 7.8), while the experimen- 
tal values aref = 71.14 A (author?) [H^l, = 64.66 A 
(author?) (Hf , and £, = 67.56 A (author?) Es- 
pecially considering the hght scattering measurements 
of the isothermal susceptibility and the turbidity of 



c. Mean and master forms of a theoretical crossover 
function 

The master modifications of the mean crossover func- 
tions of Eqs. II]) and |[2]) use the following values of three 
master (i.e., constant) factors (author?) [17] 



eii/} = 4.288 X 10^3 

L<i-^> = 25.6936 

^'{i/} = 1.73847 X 10-* 



(A8) 



The two scale factors d and ipp [see Eq. (|40|) ] needed by 
the massive renormalization scheme are related to Yc and 
Zc through the equations 



?9 = YcQ^^f'^ = 0.021069 
V'p= (Zc)"^^^i^> = 3.2507 x 10-'* 



(A9) 



In addition, the values of Eq. l(6T|) for the metric prefac- 
tors Lg £ and Xq ^ are calculated by using the following 
equations 



Z±Z± {Yc)-" 

[L{i/}x(e{i/})"] 
Z±Z± {Zc)-' {Yc)-^ 

{Z,)-\Y,)-< 

[(L{i/})"''(*{i/})"'(e{i/})^] 



0.443526 



= 0.214493 



(AlO) 

These two independent prefactors are two characteristic 
parameters of xenon, which permit to calculate all the 
other xenon prefactors of the modified crossover func- 
tions, in conformity with the two-scale-factor universal- 
ity. 



Theoretical representation of the isothermal 
compressibility in the temperature range 

r - < 5 - 10 K 



We can immediately compare all the published 
kt (AT) — i (f^) data for xenon at p — pc, to the esti- 
mated values by using the mean crossover function of Eq. 
(121) and the xenon parameter set of Eq. (|39)) . A complete 
view of the results is shown in Fig. [5] which covers about 
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Figure 5: Isothermal compressibility of xenon along the critical isochore; Only the points corresponding to the pVT data given 
in Table [V] are reported; Inserted labels are common for all the figures (see text). 



six decades on the temperature distance (the data sources 
are given belovi^; see Table rv)) . The related numerical res- 
olution of this picture is such that the size of each data 
point (including the calibration data points defined be- 
low) corresponds to a relative error-bar of ±5%, while the 
thickness of the theoretical curve (label MR) accounts 
for one of ±3%. As expected from our previous analysis 
(author?) [l3|, no difference larger than ±5% is visible 
between the curve and the experimental data points on 
the extended asymptotic domain (labeled EAD) such as 

T - < 5K, i.e., At* < L^l^ = ~ ^ x lO'^ 

(see also below the Fig. [7] and the related discussion). 



In addition, the thickness of each (colored) line, having 
characteristic slope 7e,pVT, amplitude ^~^pyrp, and color 
indexation given in Table IVj represents ±1.5% relative 
error-bar. 



a. Data sources for xenon isothermal compressibility 

As shown in the lower part of Figure \E[ the complete 
temperature range 1 mK < AT < 400 K (i.e. 3 x lO'^ < 
^'T* ^ 1) is covered by successive overlap between kt 
data which originate from three distinct experimental 
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T-Tc 


283.41 


233.41 


183.41 


133.41 


108.41 


83.41 


58.41 


33.41 


23.41 


13.41 


10.41 


8.41 


6.41 


Ref. 




0.01883 


0.02348 


0.02912 


0.04078 


0.05095 


0.06772 


0.10074 


0.18695 








0.8961 




(author?) f44] 










0.04145 


0.05185 


0.06904 


0.10227 


0.19268 


0.29065 


0.56344 




1.0099 




(author?) [47] 






















0.54791 


0.73594 




1.3020 


(author?) [12] 






















0.555 
(±0.008) 








(author?) [12] 


T~Tc 


4.91 


3.41 


2.41 


1.80 


1.4 


1.0 


0.6677* 


0.8 


0.6 


0.4 


0.2 


0.1 


0.001 












5.7371 


7.8543 


11.858 






22.701 


36.640 


79.11 


204.86 




(author?) [46] 




1.7835 


2.7224 


4.1139 










15.432 












(author?) [12] 














11.95 
(±0.15) 


19.64* 

(±1%) 












58562 
(±800) 


(author?) [12] 



Table V: Values of the isothermal compressibility kt (expressed in Pa^^) of xenon obtained from the pVT measurements along 
the critical isochore pc = 1113 kg m"^, as a function of T — Tc (expressed in K); The lower part corresponds to data obtained 
within the Ising-like extended asymptotic domain; Asterisk indicates the reference values used by Giittinger and Cannell to 
calibrate their light scattering measurements of the isothermal compressibility; (see text). 



methods: pVT measurements with labels pVT (B) (au- 
thor?) O], (M) (author?) 0, and (HS) (author?) 
[4^ : light diffusion or turbidity measurements with label 
LS (GC) (author?) [2l|, and FranhauflFer interferometry 
measurements with label IF (HM) (author?) [H,!!!]. 

Within the main part of the temperature distance to Tc 
[typically O.IK < AT < 400 K, (i.e., 3 x lO'^ < At* < 
1)], kt is obtained from pVT measurements (author?) 
[43. l4^. I47I ■ generally using graphical (author?) [iJ] or 
numerical fitting (author?) [42] of the isotherms 

to obtain the slope (^f^) at the selected critical density. 

We have also reported in Table fVl the KT-values obtained 
from pVT measurements which were used in our analysis. 
A noticeable specific situation occurs for xenon where the 
high precision of the pVT measurements of Habgood and 
Schneider (author?) ji^, can be used to check carefully 
the overlap [in the temperature range 0.2 K < AT < 
1.8 K, (i.e. 6 X 10--* < At* < 6 X IQ-^)] with the optical 
methods. 

In the intermediate temperature range [typically 
30 mK < AT < 15 K, (i.e., 10"'' < At* < 5 x lO'^)], 
the analyses of the fight (diffusion or transmission) in- 
tensity measurements (author?) j2l|, as a function of 
AT = T — Tc, provide interrelated data of the correla- 
tion length ^ and isothermal susceptibility Xp,t- kt is 

related to Xp,t by Xp.t = (^)^ = P^'^t- 

In the closest temperature range [1 mK < AT < 10 K, 
(i.e., 10"'' < At* < 3 X IQ-^)], the fiuid density profile 
versus the cell height subjected to the gravitational field 
generates Franhauffer interferograms (author?) [Hl-lsH] 
which are related to the isothermal compressibility (au- 
thor?) (ssl . [5^ . However, as already noted by Levelt- 
Sengers et al (author?) [43], the pubfished results in 
Ref. (author?) JH^I needs to be reconsidered before a 
quantitative comparison with other ones in overlaping 
temperature ranges. 

In order to focus our attention in the restricted tem- 



perature range 0.3 K < AT < 5K, around the central 
value T ~ Tc + 1.25 K (see below next §), we have also 
used the isothermal compressibility data obtained from 
the dynamic ligth scattering data of Cannell and Benedek 
(author?) [51!], and from the static ligth scattering data 
of Smith et al (author?) [H^l, as complementary data 
sources. 

A realistic estimation of the experimental uncertainty 
is of the order of 10% when comparison is made between 
the data obtained from the different methods (the optical 
method being of high relative precision but lower absolute 
precision; see for example Fig. [6] below). 

b. Effective fitting results 

To complete our analysis of the Giittinger and Can- 
nell's measurements, we have also made a comparison 
with the published fitting results, obtained from pVT 
data, using the following effective power law with an ad- 
justable non-Ising exponent, 

K*^^rl^yr{AT*r-''''^ (All) 

The je,pVT and ^tpVT values are reported in Table IVTl in 
addition with the related finite extension At* ■ — At* „^ 
of the experimental temperature range. The labels 
#6, 5, 4, of the fines of Table |Vl] correspond to the re- 
spectives labels #6, 5, 4, of the columns of Table ITll The 
typical error-bar values of these fitting results are given in 
line #6. At each restricted temperature range, we have 
attached the central value (AT*pyy) = -y/AT^^Ar^^ 
(in log scale). So that, in the next Fig. \7\ each effective 
power law result is ifiustrated as a form of a (colored) 
point-segment (with 7e,pyT-value as a label) covering the 
range At^-^^ < At* < Atj*^^, around a central (full) 
point fixed at (AT*yj,) . 

In this Table Ivj, results (author?) [l3] given in fine 
labeled #7mf correspond to a fit of the pVT data of 
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7e,pVT 


r+ 

^ e.pVT 






At* 
^ ' max 


/At* \ 


A-7-* 

(Table nil) 


-n-Zo (1 e,pVT)r+ 

E.th 


Vc 

Eq. I[74ll 


#6 


1.211 (±0.02) 


0.0743 (±0.015) 


(author?) [43] 


6.9 X lO""* 


6.2 X 10"^ 


2.07 X 10"^ 


2.95 X 10^^ 


2.3 


0.04676 


#5 


1.190 


0.0793 


(author?) [42] 


2.8 X 10"^ 


4.6 X 10"^ 


1.13 X 10"^ 


1.135 X 10"^ 


-1.9 


0.01429 


#4 


1.1665 


0.089 


(author?) [14] 


6.2 X 10~^ 


8.1 X 10"^ 


2.24 X 10"^ 


3.338 X 10"^ 


0.1 


0.02618 


#7MF 


1 


0.11 


(author?) [12] 


0.5 


1 


7.1 X 10"^ 


(oo) 


30. 


n.a. 



Table VI: Lines with labels #6, 5, 4 (corresponding to the columns with same labels in Table [11)1 : Effective exponent 7e,pVT 
(column 2) and amplitude rjf^^j, (column 3) of a power law description of Eq. l|Allll . for Aty obtained from xenon "pVT 

measurements in the temperature range and (geometrical) mean temperature (At*^j-) = y^Ar^^Ar^^ (column 6); see 
references given in column 4; Calculated values of the local temperature distance C^t^^ (column 7) are from Table [ll] linelO; 
Residuals R% (r^p^j,)^,^ (column 8), expressed in %, between experimental (r^^^j,) xenon amplitude and calculated ( F^^j^ 

) xenon amplitude from Table [III line 12; Calculated values of xenon crossover parameter "Qc (column 10) from Eq. I|74p : Line 
with label i^'yuF ■ Equivalent results for 7^ = 7mf ~ 1 (see text); n.a.: non available. 



Beattie et al (author?) [4J|, where the mean-field value 
7e = 7mf = 1 of the effective exponent can be observed at 
large temperature distance from Tc (typically At* > 0.5, 
see also Ref. (author?) [H). 

A similar introduction of the effective power law to 
analyse the results of light scattering experiments gives 
access to a quantitative comparison between the pVT 
and light scattering measurements of the isothermal com- 
pressibility. As a matter of fact, Giittinger and Cannell 
have claimed that the correction to scaHng terms are im- 
portant by demonstrating that the susceptibility deviates 
systematically from a simple power law behavior with the 
effective exponent value 7e = 1.206. More precisely. Fig. 
2 of Ref. (author?) ^ shows that 7e ~ 1.206 is the 
slope of the tangent Hne to the rough experimental be- 
havior at Ar*^^]^ 206 — 4-5 x 10"^, i.e., the temperature 
distance which corresponds to the minimum of the de- 
viation curve in this Fig. 2. We have used this result 
to renew a fit of twelve compressibility data measured in 
the restricted temperature range 9.115 x 10^* < At* < 
1.95 X 10-2 (i.e., 0.26 K < AT < 6.90 K), i.e., a tempera- 
ture range with a central value (At*qq) — 4.215 x 10~^, 
very close to the Ar*^^^ 206 o^i^. As expected, our fit- 
ting values 7e = 1.205879 and r+ = 0.07551466 (re- 
ported in line #GCe of Table lVIip . are in excellent agree- 
ment with the Giittinger and Cannell's ones 7e = 1.206 
and r+ = 0.075135. The latter amplitude value was 
calculated applying a —4% correction (corresponding to 
the fit deviation observed in this restricted tempera- 
ture range) to the effective amplitude of the published 
original fit. This resulting power law behavior may be 
seen from Fig. [6] which gives the % residual between 
the dimensionless isothermal compressibility and the fit- 
ting equation k^ ^c, = 0.07551466 (At*)"^-^''^*^^ For 
us, the most important consequence of the high rela- 
tive precision of the Giittinger and Cannell's measure- 
ments is the demonstration that a well-defined local value 
7e ~ 1.20588 of the effective exponent can be mea- 
sured at a well-defined local value {T — T^. ~ 1.260 if) 
of the temperature distance to T^.. Accordingly, by ex- 
pressing now the condition 7e,th = 1-205879 in comple- 



ment to the similar eight conditions given in Table [ill for 
the theoretical crossover function of Eq. ([2]), such a lo- 
cal value of the effective exponent must be observed at 
Ar^ {"fe.th = 1.2058) = 4.347x10-^, while the calculated 
local values of the effective theoretical and xenon ampli- 
tudes are = 0.396926 and T+^,^ {'je,th = 1.2058) = 
0.0747481. We will exploit the noticeable agreement be- 
tween the experimental and theoretical local value of the 
effective power law behavior in the last § of the Ap- 
pendix. We note that the calculated value of the cor- 



relation length is then 



;(7e~l. 20588) 



6.764, i.e, within 



our expected range for the Ising-like extended asymptotic 
domain previously defined by the condition ^* = ^ ^ 3. 

Previously, in a first approach independent of any one- 
parameter crossover theory, it is necessary to control that 
the compressibility data obtained by any other measure- 
ment method covering a similar restricted temperature 
range, are satisfying this effective power law behavior 
with comparable values of the related effective ampli- 
tudes. This is the object of the results given in lines 
#pVT, CB, SGB of Table Eni (see also Fig. [e]), where 
we have used: 

i) line i^pVT, the pVT data of Table |V] covering the 
temperature range 0.1 K < T-T^ < 10.41 K (i.e., 3.451 x 
lO"'' < At* < 3.593 x 10"^); 

ii) lines 4f^CB, the dynamic li gth scattering data of 
Cannell and Benedek (author?) j5l| covering the tem- 
perature range 0.2 K < T - < 20 K, i.e., 6.9 x lO""* < 
At* < 6.9 X 10-2, and 

iii) lines ^SGB, the static ligth scattering data of 
Smith et al (author?) [50] covering the temperature 
range 0.045 K < T - Tc < 5.1 K, i.e., 1.553 x 10"* < 
At* < 1.76 x 10-2 . 

In the pVT case, the fitting values of the exponent- 
amplitude pair confirm the power law behavior observed 
from Giittinger and Cannell's measurements of high rel- 
ative precision. On the other hand, we have reported the 
initial fitting^^ results [lines #Ci?, SGB, (a)] of Refs. (au- 
thor?) [50, 51.]. They are only used to illustrate quanti- 
tatively the effects due to a large uncertainty on measure- 
ments of the geometrical factors in light scattering exper- 
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7e 


r+ 


Ref. 


Ar* 

^ ' max 


(Ar* ) 




r+ 

e,th 


e,th 




■&C 
Eq. dZH) 


GCe 


1.205879 


0.07551466 


t.w. 


9.115 xlO"'' 
1.95 xlO"^ 


4.215 X 10~^ 


4.347 X 10"^ 


0.0747481 


1.03 


0.396926 


0.0285147 


pVT 


1.204519 


0.07590869 


t.w. 


3.451 xlO"'' 
3.593 xlO"^ 


3.521 X 10"^ 


4.778 X 10"^ 


0.075298 


0.81 


0.40195 


0.0265262 
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3.451 xlO"* 
2.5886 xlO"^ 
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(c) 


1.201586 


0.0768207 


t.w. 






5.804 X 10"^ 


0.076466 


0.46 


0.41283 


0.0238014 




L.AL ZIZ U.UO 
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1.553 xlO"'' 
1.76 xlO-2 
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(b) 


1.204381 


0.0679096 


t.w. 


1.5877 xlO-* 
1.723 xlO"^ 


1.653 X 10"^ 












(c) 


1.204381 


0.0755227 


t.w. 






4.823 X 10"-' 


0.075354 


0.22 


0.40246 


0.0224563 



Table VII: Fitting results using an effective power law equation = (Ar*) in a restricted temperature range Arj^j^, < 
Ar* < Ar,*ax (see text); The calculated local values of the temperature distance Ar,*), (column 7), the xenon amplitude F^f^j^ 
(column 8), the universal amplitude g (column 10), and the xenon crossover parameter i?£ (column 11), are for each condition 
7e,tii = 7e ; t.w.: this work. 



iments (author?) [s^l, or to an indirect estimation of the 
magnitude of the compressibiUty from an elaborate and 
complex analysis of the Brillouin spectra of xenon (au- 
thor?) j5l|. For exaraple, the numerical values found 
from Ref. (author?) [50i| with 7e = 1.21 can be multi- 
plied by a factor ]^ = ore^-'lf - 1-14 to match those 
reported in Ref. (author?) [51]. However, in spite of 
the importance (~ (10 — 20) %) of these effects, we have 
confirmed with our fitting results [lines #Ci?, SGB, (b)] 
that the value of the effective exponent is well in the range 
7e ~ 1.20 — 1.21 for the restricted temperature range se- 
lected here. In addition, the large uncertainty on the 
experimental value of the effective amplitude can be de- 
creased by using our calibrated value of the isothermal 
compressibility at T = Tc + IK. For example, after cal- 
ibration [lines #CB, SGB, (c)], the effective amplitude 
for Cannell and Benedek's data was lowered by 1.573%, 
while the one for Smith et al's data was increased by 
11,21%, in agreement with around ~ 14% initial devia- 
tion between these two data set. 

Finally, as reported in column 9 of Table IVIIj the per- 
cent deviation %i?(r+)p+ = 100 | , - 1 ) 

between the experimental and theoretical effective am- 
plitudes for the four experimental values of the effective 
exponent 7e is on the 1%-level. In Fig. [6l we have also 

reported the residuals %R {k^)qq = 100 ^ ij 

(expressed in %) between each experimental data and 
the calculated one k'^ qq^ using the effective power law 
K^p — 0.07551466 (At*) ^-^^^^'^^ ^s reference (see line 
#GCe in Table Ivnl. We note the data agreement at the 



same percent level that the precision on the calibrated 
values of the isothermal compressibility, which leads to a 
conclusion that the "experimental" values of the effective 
amplitude r+ are estimated with an uncertainty of 1%. 

Therefore, in a second approach focussed on the valid- 
ity test of the one-parameter crossover modelling in pure 
fluids, which will be discussed below in § A. 3, all these 
fitting results with effective values of the exponents ob- 
served in a small restricted temperature range at large 
temperature distance to Tc are appropriate: 

i) to check the calibration of the leading asymptotic 
amplitude F"*" (with 7 fixed); 

ii) to verify the uniqueness of the scale factor -& by 
estimating the related local values of and then, 

iii) to control the master values of the leading ampli- 
tude Z+ = 0.119 [related to r+, see Eq. ^Mil] and the 
confiuent amplitude Z^'^ = 0.555, [related to see Eq. 

mi 



c. Data comparison from reference to the master crossover 
function 

Data comparisons are magnified in the three parts (a), 
(b) , and (c) , of Figure[7] (in lin-log scale) , either as a func- 
tion of At* in the lower horizontal axis, or as a function 
of the theoretical (£th) and master {£* = ^) correlation 
lengths (author?) in the upper horizontal axis. The 
extensions of the Ising-like preasymptotic and extended 
asymptotic domains, the selected 7e-conditions of Table 
ITTl and each experimental temperature range, are illus- 
trated as in Fig. [2l 
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Figure 6: Residuals %R{k^)q^ (expressed in %) of 
the data from the effective power law fi:T,GCe ~ 

0.07551466 (Ar*)"^'^"''**''^ (see text); The horizontal and ver- 
tical (pink) lines permit to define the local values of the 
effective amplitude (r+^^ = 0.99 x 0.07551466) and the 
temperature distance [Ar/Jj (7^) = 4.347 x 10~^, see Table 
IVIII ] where the line of slope 7^ — 1.205879 and equation 
Kt = '^t,jc {At*)~''' , is tangent to the experimental singular 
behavior measured by Giittinger and Cannell, at the contact 
point of coordinates {Arj\ (70) ; 0.99k^ (J(j^ [Atj*^ (7e)]} (see 
also Table IVlTl and text); Horizontal segment labeled GCe: 
restricted temperature range of the fit by the effective power 
law KT.GCe = 0.07551466 (Ar*)"^-^"^**''^ (colored) symbols: 
see inserted legend and Tables |V] and IVIII All the other in- 
serted labels are identical to the ones previously used. 



The asymptotic behavior of the quantity 



(Ar*)- 



timated from the crossover function [(black) curve MR] 
is given in Figure [TK- Such a presentation magnifies the 
role of the theoretical 7-value on the determination of 
the leading amplitude F+close to Tc (given by an hor- 
izontal asymptotic limit for Ar* 0). Correlatively, 
at large values of Ar*, the increasing contribution of 
the confluent corrections is well demonstrated by the de- 
creasing 7e,pyT-values, here given by the direction of each 
point-segment closely tangent to the MR curve (see Ta- 
ble |Vl|. Only the published experimental data points 
[(pink) small circles] and fltting result [GC4 (pink) line] 
of Giittinger and Cannell's measurements are reported in 
this Fig. Efa). 



100 



1 



expressed 



The residual R% {k^)^^ 

in %, for all the selected data, are reported in part 
(b). The accurate (~ ±3%) theoretical representation 
of pVT data of Habgood and Schneider and light scat- 
tering data of Giittinger and Cannell confirms that the 
extended asymptotic domain is well such that Ar* < 
'^EAD '^■02 — 0.03, which well corresponds to the range 
£* = -|- > 2 — 3 of accurate theoretical representation of 
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Figure 7: Fitting results of xenon isothermal compressibility 
a function of Ar* (lower horizontal axis) or correla- 
tion lengths (upper horizontal axes), by using Eq. I|23p with 
xenon parameters of Eq. lliOl) : (a) "J^-t (full black line: 
calculated confiuent correction from the crossover function); 
(b) Residuals R%{k'^)j^^ (expressed in %), for all the ex- 
perimental data and fitting results selected in this work from 
reference to the theoretical estimation by the crossover func- 
tion; (c) Same as Fig. IDJb), calculated [line MR: using the 
scale dilatation method; line 1: from Eq. ifTTjl : line 2: from 
Eq. I|74p : colored point-segments: from Table IVll column 10; 
crosses: from Table IvTH column 11]; All the inserted labels 
and symbols are identical to the ones of the previous figures 
(see text). 
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the master singular behavior of the correlation length, as 
previously observed. 

In Figs ni&) and (b) , a significant difference occurs in 
the case of the actual mean field like behavior represented 
by the segment jpVT = 7mf = 1 (see line 5 of Table 
I]). In the temperature range 90K < AT < 300 K, (i.e., 
3 X 10~^ < At* < 1), the effective classical-to-critical 
crossover for xenon is not accounted for by the mean 
crossover function. In addition, each part (a) and (b) also 
contains the relative position of a grey segment (labeled 
vdW) which corresponds to the mean-field-like equation 
k^At* = i, where xenon is assimilated to a van der 

Waals fiuid (i.e. with 7e,vdW = 1 and ^t.vdw = 
That clearly illustrates the failure of the van der Waals 
equation of state close to the liquid-gas critical point. 

Figure [Tjc) is similar to Figure [DJb) . Now we have 
added the point-segment representation of the values of 
^CpVT calculated using Eq. I|74p and the fitting results 
of the pVT measurements with Eq. (|All|) , as reported in 
Table [Vl We will return below to the discussion of these 
pVT results which confirm that the xenon crossover is 
characterized by a unique value [Eq. (|6T|) ] of the scale 
factor d (Xe) over the temperature range where — > 
2-3. 



3. Data calibration and related discussion to the 
uniqueness of the crossover parameter 

As already evidenced in the seventy's (author?) [43 |. 
fitting the experimental singular behavior using Wegner- 
like expansion with a limited number of terms generates 
large uncertainty in the amplitude determination, due 
to the low convergence of this expansion. Examining 
more carefully the isothermal compressibility data ob- 
tained from interferometry (author?) |55,], light scat- 
tering (author?) 0, and p^T (author?) [ii, H,!^ 
measurements, it was noted by one of us (author?) |12l |. 
that a restricted three-term Wegner-like expansion, with 
fixed exponent values 7/ = 1.24±0.01 and A/ — ^±0.05, 
can provide satisfactory representation of the kt singular 
behavior in the range T — Tc < 13 K (see for example the 
curve G3 in Figure Wp) which corresponds to Eq. i(43|) 



with r+ = 0.0577, at 



1.55, 



2X 



-1.3, and 



0; see Table HVl last line). It was conjointly shown 
that r+, a+, and a^^ can also be calculated using three 
calibrated values at three selected temperature distances 
(in logarithmic scale, see below). The correlative main 
results were the estimation of min (ajmin — 0-6) and 
max (ajmax — 1-^) values of and the evaluation of 
the correlated error-bars which can then be controlled 
by a careful analysis of the residuals. 

We recall the main advantages of such a three points 
calibration approach, first by examining the estimation 
of the leading amplitude and its attached uncertainty for 
the closest point to the critical point. 



Calibration of the leading amplitude the closest to the 
critical point 



At T - Tc 



ImK (Ar^ 



ci 



3.45137 X 10-6), 



which corresponds to the "lowest" temperature distance 
accessible by interferometry experiments of Hocken 
and Moldover (HM) (author?) fH^, the standard 
dimensionless value of the isothermal compressibility 



was defined 



r+ 

^ MR6 



(Ar^i) 



such that 

7MR6m, 



'^T.l (^'''Cl) 



(Ar5i)' 



10^ 
and 



were 
with 



= (3.415 ±0.035) X 
with {7/ = 1.240; r+ = 0.0574 ± 0.0006} 
{7MR6„„ = 1.24194; r+j,g_ = 0.0563 ± 0.0006}. 
These initial values of the leading amplitude r+ 
re-evaluated accounting for the small difference 
7 = 1.2395935. Similar re-evaluation was conjointly 
made for the leading amplitudes (of the dimensionless 
correlation length ^*) and B (of the dimensionless order 
parameter density Ap^y) due to the small differences in 
V and l3 values. The updated values of the corresponding 
leading amplitudes (such as r+ (Xe) = 0.05782 ± 0.0006) 
were then used to optimize the new central values of the 
corresponding master amplitudes (such as Z+ — 0.119 
reported here), using Eqs. (|A6p . corresponding values 
of the critical parameters [see Eq. (jA2p ]. and the 
updated value of the universal amplitude combina- 



tion 



Z+ 



-d 



R+ = 2.92638 (with R+ = 0.2696967 and 



R 



c 
Z+ 



-- 0.057406 (author?) [7]). Finally, selecting 
^ - 0.119 and Zm = 0.468 (author?) Ill as "central" 
values for two independent master amplitudes, any 
other (master and physical) leading amplitude can be 
calculated with the numerical precison of the massive 
renormalization scheme. For example, the number of 
digits in the quoted value = 0.570481 is similar to 



the numerical 



precision of ZJ = 2.121008 (see Table 
IBi), and more generally, data of Eqs. (|A6P to (|Alip are 
obtained according to this scheme. 

Considering now the 7 differences between the theo- 
retical estimations, we can admit that the mean value 
7 = 1.2395935 (see Table [lb) is on the range 710W — 
1.235 < 7 < 7high ^ 1.242. The max value 7high 1-242 
was used in our initial fitting of the GC data (see line 
2 of Table [ml, while the min value 710W — 1.235 was 
recently estimated from the minimal subtraction scheme 
(author?) [2§|. The correlative effect on the r"*" value is 
of the order of ±3%, as shown by the curves (labeled 
7high and 710W in Fig. [Tj:) related to the power law 



'^high.low 

and r+„ 



high,iow(Ar*)-^-^°"", with FjTigh = 0.056 
0.062, respectively. This theoretical uncer- 
tainty level appears then comparable to the experimen- 
tal one (~±5%), when the 7 and values are obtained 
from interferometry measurements close to the critical 
point. For example, the initial "free" values are 7 = 1.23 
and r+ = 0.062 ± 0.006 (author?) fH^I, while a subse- 
quent analysis made for fixed 7 = 1.24 gives r+ = 0.058± 
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0.002 (author?) [s^. As a result, we note that the cal- 
ibrated value of 1 (Ar^i) = (3.415 ± 0.035) x 10^ re- 
mains well representative of the theoretical analyses of 
the interferometry measurements at this closest temper- 
ature distance to Tc. However, we also recall that the 
uncertainty associated to the experimental determination 
of the critical temperature (roughly estimated of the or- 
der of ±0.5 mK) cannot be accounted for in the above 
analyses of the interferometry measurements. 

To partly conclude, when the exponent 7 is fixed 
to its Ising-like theoretical value 7 = 1.2395935, the 
above standard value of the isothermal compressibility 
at r — Tc = 1 mK is a reahstic asymptotic constraint to 
fix the "central" value r+ (Xe) = 0.05782(±0.0006) of the 
leading amplitude, before to analyse the isothermal com- 
pressibility data over larger temperature distances to the 
(known) critical temperature with the objective to es- 
timate the contribution of the confiuent corrections to 
scaling in the xenon case. 



gives a+ = af^ = 1.25557, while Eq. IjAlSP gives 
aj^ = —0.569934. The small difference with the value 
a J — 1.22961 calculated using Eq. l(39l). is due to the 
similar "constrained" adjustment of the first amplitude 
a'lj of the order parameter density (not reported here) , 

which maintains the universal value -¥- = 0.9 of the re- 

lated amplitude ratio (author?) [9]. In such an opti- 
mized result from two properties, the best central values 
of the master confiuent amplitudes are -Z^'"*" = 0.555 and 
~ 0-4995, respectively. However, as previously indi- 
cated, a true error-bar of the order of ±35% (at least) can 
be attached to these confiuent amplitudes in the absence 
of data calibration closer to Tc. 



c. The uniqueness of the crossover parameter and the 
effective extension of the critical domain beyond the 
Ising-like preasymptotic domain 



b. Calibration of the first amplitude of the critical confluent 
correction to scaling beyond the Ising-like preasymptotic 
domain 

As indicated above, by using two other points of cali- 
bration properly selected to cover the temperature range 
of optical and pVT measurements, provides an ana- 



lytical estimation of a+ 



and aj^ in Eq. 



[with 7, A fixed, and a^^ = 0]. Indeed, from the cali- 
brated values { At(^2 ; 2 } ^-n^ { At^s ; '^t 3 } ) can de- 

1 y and 



fine the pairs Atj 

ArS3;il3 = 



R32 ' 2 



r+(Ara3)"" 



1 



r+(Ar£,)- 

We obtain 



Y2 {R32f^ - Y3 



Ar, 



C3 



(i?32)' 



(A12) 



and 



Y3 - Y2 (i?3 



n2A 



*2x 



(Ara^Ar^g) 



(i?32)^ - 1 



(A13) 



Equations (|Al2p and I|A13P have convenient analytic 
forms to check the infiuence of the selected values for 
the exponents and calibrated points (author?) [l^ . 
For example, we can choose the caHbrated points at 
T2 - Tc = IK (Ar^j = 3.45137 x 10"^) with kt,2 = 
(1.195 ±0.012) X 10~5pa"^ {k*^ ^ = 69.8 ± 0.7), and 
Ts-Tc^ 13.41 K (Ar^g = 4.62829 x 10"^) with kt,3 = 
(5.55 ±0.08) X 10-^Pa~^ (^^ 3 = 3.24 ±0.04) (see Table 
rv|) . Using the updated theoretical values 7 = 1.2395935, 
A = 0.50189 (see Table [lb), r+ (Xe) = 0.0578204, and 
the two calibrated pairs {i?32 = 13.41; Fa = 0.071051} 
and { At53 = 4.62829 lO'^; Y3 = 0.242478}, Eq. ([Anil 



Since the asymptotic amplitude r+ is now fixed, the ef- 
fective exponent-amplitude pair 7e,pV'T; ^tpVT reported 
in columns 2 and 3 of Table IVI[ can be used to calcu- 
late the value of a crossover parameter dc,pVT at each 
"local" temperature distance (AT*y^), using the scale 
transformation of Eq. (l73l) . and the theoretical values 



of Z+ and Z+g (7e,pyT) given in Tables [J and [TTl The 



corresponding values of 'dc,pVT\ 



lpVT=cte 



are given in last 



column of Table [VB Each pVT result for 'dc,pVT\^ 
as a function of At* , has then been illustrated in Fig. [3 
(c) as a point-segment form. The "local" value (Ar*yj,) 
is close to the calculated one Ar^*;^ (see columns 6 and 
7 in Table [Vj). In addition, the %-residuals = 

^ c,th 

100 I — 1 I between the effective amplitudes 

compare favorably with the estimated experimental pre- 
cision (see column 8 in Table [Vl[) . Therefore, at large 
temperature distance from Tc, we obtain a significant 
confirmation that the values of the effective crossover pa- 
rameter are close to the one of the asymptotic scale factor 
d = 0.021069, and can then be considered as independent 
of Ar* in the Ising-Hke extended asymptotic domain. 

Using in a similar manner the fitting results given 
in Table IVIIi we have obtained the corresponding val- 
ues of dc given in column 11. However, we have now 
a better control of the related uncertainty, thanks to 
the high precision of the ligth scattering experiment of 
Giittinger and Cannell. As a matter of fact, admit- 
ting in a first approach that the values of F"*", (Ar*^), 
and Z+g (7e) are known with zero uncertainty, while 

= rc+th'(i + ^r),with 5t ^ 

from Eq. ((74|) and d = 



1 — 2%, it is easy to show 



r+ 



that 



{i + 5Ty 



ST 



(A14) 
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Looking then to the percent deviation reported in Fig. 
El we can observe that the true tangent (pink) Hne 
of slope 7e = 1.20588 has an ampUtude ~ 1% lower 
than the amplitude of the effective power law — 

0. 07551466 (Ar*)-'-^^^*'^ used as a reference (see also 
Table IVIII column 9, where the calculated value of the 
residual is -1% in line #GCe) . From Eq. I|A14P with 
5T = 0.01 and dc ^ 0.02851 (see line #GCe, colum 
11, Table En]), we obtain i) ~ 0.021988 which is in ex- 
cellent agreement (+4.38%) with our initial estimation 
1? ~ 0.21069 from the scale dilatation method. As previ- 
ously underlined, the precise description by a local expo- 
nent value defining the slope of the tangent line to the sin- 
gular behavior of the isothermal compressibility of xenon 
at a well-defined temperature distance to Tc, is one of the 
major points of interest of the Giittinger and Cannell's 
results to validate the one-parameter crossover modelling 
predicted by the massive renormalization scheme. In 
a similar manner, using Eq. I|A14P with ST and 
given in lines #pVT, CB(c), SGB(c) of Table Elll we 
obtain 1? ~ 0.02160, 0.021232, 0.021136, and the corre- 
sponding deviations-l-2.52%, 0.77%, 0.32%), respectively 
for the three other fitting results reported in Table IVIII 

Obviously, similar effect is produced by the uncertainty 
level attached to the determination of the leading ampli- 
tude r+, justifying oncemore its independent estimation 
from a "standard" value of the isothermal compressibil- 
ity at a temperature distance well-inside the Ising-like 
preasymptotic domain. 

Henceforth, the importance of the scaling form of Eq. 
([74l) in the determination of a unique asymptotic value 
for the scale factor ■& is clearly established. Moreover, 
Eq. ([Tl]) is vaHd in the range ^ > 2.5 - 3, or At* < 
£g^Q ~ 0.02 in xenon case. 

One complementary remark can be formulated. 

Outside the Ising-like extended asymptotic domain, 

1. e., typically for At* > 10^^ in xenon case, it is well- 
established that the Ising-like universality is not vaHd. 



For example, increasing the temperature distance to Tc 
on the pVT data analyzes, we are able to observe the con- 
tinuous decreasing behavior of je,pVT until a value close 
to mean-field value je,pVT ~ 7mf — 1 when At* ^ 1 (see 
Table IVip . That unambiguously discriminates a sharp 
domain, i.e. typically 0.3 < At* < 0.5, where 7e,pyT 
crosses the "mean" crossover value 71 = '^'^^'^^ « 1.12 
(as initially reported in Ref. (author?) (2^). For xenon, 
this classical-to-critical crossover "crossing" is expected 
close to < 1, that means that the correlation length is 
of the same order of magnitude or lower than the short- 
range molecular interaction. In Table lVTl the non-defined 
value i}c,pVt\^ vt=i ^'^ corroborates that the mean- 
field behavior predicted by the theoretical crossover func- 
tion with 1? = 0.0210, is not compatible with the pVT 
experimental result F+^j^^ ~ 0.11. Such a typical limit 
of the mean-field-Hke range (see upper part of Figure 
[si is well-illustrated by the corresponding transforma- 
tion between the opened squares represented in Figure [21 
which cannot account for the "experimental" location of 
the opened circle when 7e = 7mf — 1- 

Finally, Figure |7| (c) confirms that the effective ex- 
tended critical domain of xenon, corresponding to the 
condition ^ ^ 3 discussed in a detailed manner in 
Ref. (author?) [l3|, is well characterized by a single 
crossover parameter whose value at the largest temper- 
ature range is comparable to the one of the asymptotic 
scale factor estimated from the scale dilatation method. 
We can conclude that singular behavior of any property 
in this extended critical domain of xenon can be calcu- 
lated in conformity with the universal features predicted 
by the massive renormaHzation scheme, only using the 
required four critical coordinates to define the position 
and tangent surface of its actual liquid-gas critical point 
on the p, Vp, T phase surface, as expected by one of us 
two decades ago. 



[1] M. A. Anisimov and J. V. Sengers, in Equations of State 

for Fluids and Fluid Mixtures, Part I, J.V. Sengers, R.F. 

Kayser, C.J. Peters, and H.J. White, Jr., Eds. (Elsevier, 

Amsterdam, UK, 2000) pp. 381-434. 
[2] see for example Zinn Justin, Euclidean Field Theory and 

Critical Phenomena , 3'^^ ed. (Oxford University Press, 

1996). 

[3] A. Pellissetto, P. Rossi, and E. Vicari, Phys. Rev. E 58, 
7146 (1998). 

[4] C. Bagnuls and C. Bervillier, J. Phys. (Paris) Lett., 45, 
L-95 (1984). 

[5] C. Bagnuls and C. Bervillier, Phys. Rev. B 32, 7209 
(1985). 

[6] C. Bagnuls and C. Bervillier, D. Meiron, and B. Nickel, 
Phys. Rev. B 35, 3585 (1987); 65, 149901(E) (2002). 
Verifier la page 3585. 

[7] C. Bagnuls and C. Bervillier, Phys. Rev. E 65, 066132 



(2002). 

[8] R. Guida and J. Zinn- Justin, J. Phys. A: Math. Gen. 31, 
8103 (1998). 

[9] Y. Garrabos and C. Bervillier, Phys. Rev. E 74, 021113 
(2006). 

[10] Y. Garrabos, F. Palencia, C. Lecoutre-Chabot, C. J. 
Erkey, and B. Le Neindre, Phys. Rev. E 73, 026125 
(2006). 

[11] Y. Garrabos, C. Lecoutre-Chabot, F. Palencia, D. 

Brosetta, B. Le Neindre, and C. J. Erkey, Phys. Rev. 

E 75, 061112 (2007). 
[12] Y. Garrabos, Ph. D. Thesis, University of Paris (1982). 
[13] Y. Garrabos, J. Phys. (Paris) 46, 281 (1985) [for an en- 

glish version see e-print/cond-mat/05 12408]. 
[14] Y. Garrabos, J. Phys. (Paris) 47, 197 (1986). 
[15] Y. Garrabos, B. Le Neindre, R. Wunenburger, C. 

Lecoutre-Chabot, and D. Beysens, Int. J. Thermophys. 



29 



23, 997 (2002). 
[16] Y. Garrabos, Phys. Rev. E 73, 056110 (2006). 
[17] Y. Garrabos, C. Lccoutrc-Chabot, F. Palcncia, B. Lc 

Neindre, and C. J. Erkey, preprint (2006) (see cond- 

mat/). 

[18] F. J. Wegner, Phys Rev. B 5, 4529 (1972). 

[19] K. G. Wilson and J. Kogut, Phys. Rep. 12 C, 75 (1974). 

[20] In this approach of scaling, we use the correlation length 
and the susceptibility as independent properties in terms 
of the two-scale-factor universality. We can then ignore 
the universal features related to the additive univer- 
sal constant of the heat capacity. Moreover, in a pure 
fluid case, the critical background constant is mixed with 
the non-universal contribution of the regular background 
terms. 

[21] H. Giittingcr and D. S. Cannell, Phys. Rev. A 24, 3188 
(1981). 

[22] C. Bagnuls, C. Bervillier, and Y. Garrabos, J. Phys. 

(Paris) Lett. 45, L-127 (1984). 
[23] C. Edwards, J. A. Lipa, and M. J. Buckingham, Phys. 

Rev. Lett. 20, 496 (1968). 
[24] J. M. H. Levelt Sengers and J. V. Sengers, in "Progress in 

Liquid Physics", Ed. C. A. Croxton (John Wiley & Sons, 

New York, 1978) pp. 103. 
[25] V. Privman, P. C. Hohenberg, and A. Aharony, Universal 

critical point amplitude relations, in "Pfease Transitions 

and Critical Phenomena", Vol. 14, Ed. C. Domb and J. 

B. Lebowitz (Academic Press, New York, 1991). 
[26] M. A. Anisimov, A. Povodyrev, V. Kulikov, and J. V. 

Sengers, Phys. Rev. Lett. 75, 3146 (1995). 
[27] E. Luijten and H. Meyer, Phys. Rev. E 62, 3257 (2000). 
[28] L Hahn, F. Zhong, M. Barmatz, R. Haussmann, and J. 

Rudnick, Phys. Rev. E 63, 055104(R) (2001). 
[29] F. Zhong, M. Barmatz, and L Hahn, Phys. Rev. E 67, 

021106 (2003). 

[30] F. Zhong and M. Barmatz, Phys. Rev. E 70, 066105 
(2004). 

[31] V. Dohm, Z. Phys. B: Condens. Matter 60, 61 (1985); 
R. Schloms and V. Dhom, Europhys. Lett. 3, 413 (1987); 
Nucl. Phys. B 328, 639 (1989); Phys. Rev. B 42, 6142 
(1990). 

[32] E. Luijten, H. W. J. Blote, and K. Binder, Phys. Rev. 

Lett. 79, 561 (1997). 
[33] E. Luijten and K. Binder, Europhys. Lett. 47, 311 (1999). 
[34] E. Luijten and K. Binder, Phys. Rev. E 58, R4060 (1998); 

59, 7254(E) (1999). 
[35] M. H. Miiser and E. Luijten, J. Chem. Phys. 116, 1621 

(2002). 

[36] Z. Y. Chen, P. C. Albright, and J. V. Sengers, Phys. Rev. 

A 41, 3161 (1990). 
[37] Z. Y. Chen, A. Abbaci, S. Tang, and J. V. Sengers, Phys. 

Rev. A 42, 4470 (1990). 
[38] M. A. Anisimov, S. B. Kiselev, J. V. Sengers, and S. 

Tang, Physica A 188, 487 (1992). 
[39] Belyakov and Kiselev, Physica A 190 (1992). 
[40] V. A. Agayan, M. A. Anisimov, and J. V. Sengers, Phys. 

Rev. E 64, 026125 (2001). 
[41] J. S. Kouvel and M. E. Fisher, Phys. Rev. 136, A 1626 

(1964). 

[42] J. M. H. Levelt Sengers and J. V. Sengers, Phys. Rev. A 

12, 2622 (1975). 
[43] J. M. H. Levelt Sengers, W. L. Greer, and J. V. Sengers, 

J. Phys. Chem. Ref. Data 5, 1 (1976). 
[44] J. A. Beattie, R. J. Barriault, and J. S. Brierley, J. Chem. 



Phys. 19, 1219 (1951); 19, 1222 (1951). 
[45] M. A. Weinberger and W. G. Schneider, Can. J. Chem. 

30, 422 (1952); 30, 847 (1952). 
[46] H. W. Habgood and W. G. Schneider, Can. J. Chem. 32, 

98 (1954); 32, 164 (1954). 

[47] A. Michels, T. Wassenaar, and P. Louwerse, Physica 20, 

99 (1954). 

[48] V. A. Rabinovich, L. A. Tokina, V. M. Bcrczin, Tcplo. 

Vys. Temp. 11, 64 (1973); V. A. Abovskii and V. A. 

Rabinovich, Teploenerg. 3, 44-57 (1971); 6, 45-93 (1973). 
[49] M. Giglio and G. B. Benedek, Phys. Rev. Lett. 23, 1145 

(1969). 

[50] L W. Smith, M. Giglio, and G. B. Benedek, Phys. Rev. 

Lett. 27, 1556 (1971). 
[51] D. S. Cannell and G. B. Benedek, Phys. Rev. Lett. 25, 

1157 (1970). 

[52] H. L. Swinney and D. L. Henry, Phys. Rev. A 8, 2566 

(1973) . 

[53] H. Giittinger and D. S. Cannell, Phys. Rev. A 22, 285 
(1980). 

[54] W. T. Estler, R. Hocken, T. Charlton, and L. R. Wilcox, 

Phys. Rev. A 12, 2118 (1975). 
[55] R. J. Hocken and M. R. Moldover, Phys. Rev. Lett. 37, 

29 (1976). 

[56] J. V. Sengers and M. R. Moldover, Phys. Lett. 66A, 44 
(1978). 

[57] W. J. Camp, D. M. Saul, .J. P. Van Dyke, and M. Wortis, 

Phys. Rev. B 14, 3990 (1976). 
[58] D. S. Gaunt and M. F. Sykes, J. Phys. A 12, L25 (1979). 
[59] G.A. Baker, B. G. Nickel, M. S. Green, and D. L Meiron, 

Phys. Rev. Lett. 36, 1351 (1976); G.A. Baker, B. G. 

Nickel, and D. L Meiron, Phys. Rev.B 17, 1365 (1978). 
[60] J. C. Le Guillou and J. Zinn- Justin, Phys. Rev. Lett. 39, 

95 (1977). 

[61] For example, the microscopic spacing order in xenon 
fluid can be driven by the Lennard Jones like diame- 
ter a^-' of xenon atom, which typically is cr^^ (Xe) ~ 
0.4 nm, i.e. comparable in order of magnitude to the 
estimated value Lq,/; (Xe) ~ 0.391 nm. However, this 
typical size of xenon atom is one order of magnitude 
greater than [go,c (Xe)]^^ ~ 34.3 10^'^ nm and also sig- 
niflcantly greater than the critical thermal wavelength 
At,c (Xe) ~ 8.95 10"^ nm. 

[62] L. M. Staccy, B. Pass, and H. Y. Carr, Phys. Rev. Lett. 
23, 1424 (1969). 

[63] A. B. Cornfeld and H. Y. Carr, Phys. Rev. Lett. 29, 28 
(1972); 29, E320 (1972). 

[64] V. G. Baidakov, A. M. Rubshtein, V. R. Pomortsev, and 
L J. Sulla, Phys. Lett. A 131, 119 (1988). 

[69] C. E. Hayes and H. Y. Carr, Phys. Rev. Lett. 39, 1558 
(1977). 

[66] C. W. Garland and J. Thoen, Phys. Rev. A 13, 1601, 

(1976) . 

[67] J. Thoen and C. W. Garland, Phys. Rev. A 10, 1311, 

(1974) . 

[68] T. S. Chang, C. W. Garland, and J. Thoen, Phys. Rev. 

A 16, 446, (1977). 
[69] C. E. Hayes and H. Y. Carr, Phys. Rev. Lett. 39, 1558 

(1977) . 

[70] D. Balzarini and O. G. Mouritsen, Phys. Rev. A 28, 3515 

(1983). 

[71] U. Narger and D. A. Balzarini, Phys. Rev. B 42, 6651 
(1990). 

[72] K. A. Gillis, L L Shinder, and M. R. Moldover, Phys. 



30 



Rev. E 70, 021201 (2004); 72, 051201 (2005). 
[73] As a typical example related to the master value Zm = 
0.468, the leading amplitude B = Zm (^c)~^ [Yc)'^ (with 
P — 0.3258) of the singular behavior of the order param- 



eter density of xenon takes the value B = 1.46998 which 
compare remarkably to the more recent experimental val- 
ues B = 1.470 ± 0.010 and B = 1.479 ± 0.011 (with 
/3 = 0.327) obtained in Ref. (author?) (nil- 



